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Abstract 

Splitting trees are those random trees where individuals give birth at constant rate during 
a lifetime with general distribution, to i.i.d. copies of themselves. The width process of 
a splitting tree is then a binary, homogeneous Crump-Mode-Jagers (CMJ) process, and 
is not Markovian unless the lifetime distribution is exponential (or a Dirac mass at {oo}). 
Here, we allow the birth rate to be infinite, that is, pairs of birth times and lifespans of 
newborns form a Poisson point process along the lifetime of their mother, with possibly 
infinite intensity measure. 

A splitting tree is a random (so-called) chronological tree. Each element of a chrono- 
logical tree is a (so-called) existence point (v, t) of some individual v (vertex) in a discrete 
tree, where r is a nonnegative real number called chronological level (time). We introduce 
a total order on existence points, called linear order, and a mapping ip from the tree into 
the real line which preserves this order. The inverse of tp is called the exploration process, 
and the projection of this inverse on chronological levels the contour process. 

For splitting trees truncated up to level r, we prove that thus defined contour process is 
a Levy process reflected below r and killed upon hitting 0. This allows to derive properties 
of (i) splitting trees: conceptual proof of Le Gall-Le Jan's theorem in the finite variation 
case, exceptional points, coalescent point process, age distribution; (ii) CMJ processes: one- 
dimensional marginals, conditionings, limit theorems, asymptotic numbers of individuals 
with infinite vs finite descendances. 
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1 Introduction 



A splitting tree [TT] is a tree formed by individuals with i.i.d. lifespans, who give birth at the 
same constant rate, while they are alive, to copies of themselves. The number of individuals 
alive at time t evolves with t according to non-Mar kovian dynamics, known as the Crump- 
Mode-Jagers (CMJ) process. Actually, general CMJ processes feature birth intensities that 
may vary through the lifetime of an individual as well as clutches of random sizes. Thus, the 
CMJ processes considered here are actually binary homogeneous CMJ processes. On the other 
hand, we consider splitting trees that are more general than those considered in [llj . since 
here individuals may have infinitely many offspring. This is done by assuming that for each 
individual, the birth times and lifespans of her offspring form a Poisson point process with 
intensity dtA(dr), where A is a Levy measure on (0, oo] called the lifespan measure, which has 
/ °°(1 Ar)A(dr) < oo. 

Contour processes of splitting trees as defined here are also different from those in [TT]. Our 
contour process (X s , s > 0) visits once and once only all the instants at which any individual of 
the tree is alive, that we call existence points (the set of existence points of a given individual 
is merely its lifetime, which is in bijection with an interval of the real line). To avoid confusion 
with the usual denomination of 'time' given to the index of stochastic processes (viz. the con- 
tour process in the present setting), we will call level, or chronological level, the real, physical 
time in which live the individuals of the tree. With the same goal, we will also try and give a 
Greek letter (such as r or a) to these levels, and a Latin letter (such as t or s) to time as index 
of stochastic processes. In the case when each individual has finitely many offspring, one can 
set a rule for the contour as follows: when the contour process is about to visit a birth level, it 
jumps to the death level of the newborn and then decreases linearly (at speed —1) along the 
lifetime of this individual until it encounters a new birth event. When the contour process ends 
up its visit of an individual, its value is thus the birth level of this individual. It then continues 
the visit of its mother's lifetime at the level where (when) it had been interrupted. Hereafter, 
this process will be called jumping chronological contour process, abbreviated as JCCP. 

The key result of the present work is that thus defined contour process X for splitting trees 
is a Levy process. Our inspiration comes from a previous study [24J in the critical case with 
exponential lifespans (the only case when the CMJ process is Markovian, except the Yule case, 
where lifespans are a.s. infinite). 

Our result yields a new interpretation of a now famous connection between Levy processes 
and branching processes [20, 9, 14J that can be worded as follows: the genealogy of (continuous- 
state) Markov branching processes can be coded by a scalar Markov process (called genealogy- 
coding process in p3j), which is a Levy process in the subcritical and critical cases. This 
process is usually an abstract object that has to be considered ex nihilo, whereas here, it is 
defined as the contour process of a predefined tree (another interpretation was given in terms 
of queues in |20j). 

This interpretation of a spectrally positive Levy process as a contour process living in the 
space of chronological levels of a splitting tree yields an elegant way of considering and inferring 
properties related to the individuals alive at a fixed level r (number, coalescence levels, ages,...) 
in the tree. 

In addition, it relates for the first time the genealogies of continuous-state branching processes 
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defined in an apparently unrelated way in [20] and [6j, in the finite variation case. In the former 
work, the starting object is a Levy process X with no negative jumps starting from \ an d killed 
when it hits 0, and to each timepoint t is associated a height H t (generation, integer number), 
which is given by the following functional of the path of X 



The Lebesgue measure of the set of timepoints with height n is a nonnegative real number Z n , 
and (Z n ;n > 0) is shown to satisfy both the Markov property and the branching property. 
Here and elsewhere, we will say that a stochastic process Z with nonnegative values satisfies 
the branching property, if for any two independent copies Z' and Z" of Z started respectively 
at x and y, Z' + Z" has the law of Z started at x + y. 

In the latter work, the starting object is a sequence of i.i.d. subordinators with zero drift, 
S\, S2, ■ ■ ■, and the population size at generation n is Z' n := S n o • • • o Si(x), with Z' Q = \. The 
genealogy of this continuous population is defined as follows: a point c in generation n (c £ 
[0, Z' n ]) is the daughter of a point b in generation n — 1 (b G [0, Z' n _-^\), if S n (b—) < c < S n (b). 
In particular, all points of generation n — 1 that are no jump times of »!?„ have no descendance. 

Actually, these two genealogies can be coupled simultaneously starting from a single random 
object, namely, a splitting tree T: the aforementioned Levy process X is actually the JCCP 
(jumping chronological contour process, as defined previously) of T, and Z n = Z' n is the sum 
of all lifespans of individuals belonging to generation n. Let us try and explain this briefly. 
First, at time t, the JCCP X visits (the existence point at real time Xt of) an individual whose 
generation in the discrete tree is the integer Ht, which can be seen to be exactly as in ([1]) (see 
forthcoming Corollary 13 . 5|> . Second, the time Z n spent by the contour X at height n, which 
is the total time spent by the height process H at n, is exactly the Lebesgue measure of the 
set of existence points, in the splitting tree, of individuals of generation n, which is exactly the 
sum of their lifespans. Third, one can embed the existence points of individuals of generation 
n into the real half-line by arranging lifetimes 'end to end' on the interval [0,Z n ]. Then by 
construction of the splitting tree, Z n+ \ is the sum of all atoms of i.i.d. Poisson measures with 
intensity A defined on lifetime intervals whose lengths sum up to Z n . As a consequence, condi- 
tional on Z n = z, we have obtained that Z n+ \ = S n+ \(z), where the (S^; i > 1) are independent 
subordinators with zero drift and Levy measure A (see also the proof of Theorem I4.ip . This 
explains why Z = Z', and that the genealogy defined thanks to these subordinators is exactly 
that associated with the topology of the splitting tree. 

In his seminal article |13j . M. Jirina introduced for the first time (multidimensional) Markov 
processes with continuous-state space satisfying the branching property. He studied the whole 
class of such processes in discrete time, but only a subclass of them in continuous time (pure- 
jump processes). For this reason, and because the term 'Jirina process' for continuous-state 
branching processes in continuous time, has progressively disappeared in the last fifteen years 
(in favor of 'CSBP', or 'CB-process'), I propose to call Jirina processes those branching pro- 
cesses in discrete time and continuous state-space like (Z n ;n > 0). 

Last, to clear up the difference between CMJ processes and Markov branching processes 
such as Jirina or Bienayme-Galton-Watson (BGW) processes, note that both count the 'pop- 
ulation size' as time runs, but the former count the number of individuals alive at the same 
chronological level, whereas the latter count those indexed by the successive generations of the 
tree. Replacing time by generations guarantees their Markov property to branching processes, 



H t := Card{0 < s < t : X s _ < inf X r }. 




3 



but makes their genealogy harder to infer [231 [T5] . 

In the next section, we set up classical notation on discrete trees and define plane chronological 
trees, that we endow with a genealogical structure, a distance, a linear order, a Lebesgue mea- 
sure, and a closure. We also recall well-known facts about spectrally positive Levy processes 
as well as Jirina processes. 

In Section 3, we introduce and study an order-preserving bijection ip between the closure of 
a chronological tree with finite Lebesgue measure, and a compact interval of the real line. The 
inverse if" 1 of this bijection is called the exploration process. It is its projection on chronological 
levels P2 o tp" 1 , which is called jumping chronological contour process, or JCCP. 

It is only in Section 4 that we consider random chronological trees, called splitting trees (as 
defined earlier). We study the properties of a splitting tree, and prove that the JCCP of its 
truncation up to level r is a Levy process reflected below r and killed upon hitting 0. 

This last result allows to derive a number of properties of splitting trees and CMJ processes, 
which is done in the last section. For splitting trees, an intuitive proof of Le Gall-Le Jan's 
theorem is given in the finite variation case (cf. above); the set T of levels where the population 
size is infinite is considered: when A is finite, T is empty, and when A is infinite, T has zero 
Lebesgue measure but is everywhere dense a.s.; the coalescence levels of individuals alive at the 
same level are shown to be independent and with the same distribution, that we specify; the law 
of ages and residual lifetimes of individuals from a same level is computed. For CMJ processes, 
the one-dimensional marginal is shown to be modified geometric; the supercritical CMJ process 
conditioned on extinction is characterized; various limit theorems are given, among which the 
convergence in distribution, conditional on non-extinction, of the numbers of individuals with 
infinite vs finite descendances to (p£, (1 — p)(,), where £ is an exponential random variable with 
parameter p = 1 — J °° re~ r,r A(dr), rj being the Malthusian parameter. 

Actually, the set of points with infinite descendance, or skeleton, has a discrete branching 
structure (that of a Yule tree). In a forthcoming work |17j . we extend the study of splitting 
trees with finite variation to splitting trees with infinite variation, and prove in particular that 
the branching structure of the skeleton is again discrete (see also [5]). Note that in the infinite 
variation case, lifespans are not even locally summable, and that in the presence of a Brownian 
component, we cannot stick to the tree space we deal with here. We also point out that an 
application of the present work to allelic partitions is available [16] . 

The reader who might like to avoid technicalities can proceed as follows: in Section 2, just 
focus on the definitions and the two statements; restrict the reading of Section 3 to Theorem 
13.11 Definition 13.21 and Theorem 13.3) skip proofs in the last two sections. 

2 Preliminaries on trees and stochastic processes 
2.1 Discrete trees 

Let N denote the set of positive integers. Locally finite rooted trees \21\ [9] can be coded thanks 
to the so-called Ulam-Harris-Neveu labelling. Each vertex of the tree is represented by a finite 
sequence of integers as follows. The root of the tree is 0, the j-th child of u = (ui, . . . , u n ) E N n , 
is uj, where vw stands for the concatenation of the sequences v and w, here uj = {u\, . . . , u n , j). 
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Set |w| = n the generation, or genealogical height, of u. More rigorously, let 

oo 

U={jN n 

n=0 

where N° = {0}. A discrete tree T is a subset of U such that 

(0 € r 

(ii) if v = uj G T, where j € N, then u G T 

(iii) for every u G T, there is a nonnegative integer < oo (the offspring number) such 
that uj G T if and only if j G {1, . . . , 

Note that individuals can have infinitely (but countably) many offspring. We write u -< u if u 
is an ancestor of u, that is, there is a sequence w such that u = uw. For any u = (ui, . . . , u n ), 
u\k denotes the ancestor (u\, . . . , Uk) of u at generation fc. We denote by u A v the most recent 
common ancestor, in short mrca, of u and u, that is, the sequence w with highest generation 
such that w -< u and uu -< v. 

2.2 Chronological trees 

Chronological trees are particular instances of M-trees. For further reading on M-trees, see e.g. 
|10| [TJ [19] and the references therein. The M-trees we consider here can roughly be seen as the 
set of edges of some discrete tree embedded in the plane, where each edge length is a lifespan. 

Specifically, each individual of the underlying discrete tree possesses a birth level a and a 
death level uj, both nonnegative real numbers such that a < to, and (possibly zero) offspring 
whose birth times are distinct from one another and belong to the interval (a, uj). We think 
of a chronological tree as the set of all so-called existence points of individuals (vertices) of the 
discrete tree. See Fig. [1] and [2] for graphical representations of a chronological tree. 

Definition. More rigorously, let 

\] = U x [0,+oo), 

and set p := (0,0). 

We let pi and p2 stand respectively for the canonical projections on U and [0, +oo). 

The first projection of any subset T of U will be denoted by T 

T := Pl (T) = {u G T : 3a > 0, (u, a) G T}. 

A chronological tree T is a subset of U such that 

(i) p G T (the root) 

(ii) T is a discrete tree (as defined in the previous subsection) 

(iii) for any u G T, there are < a(u) < oj(u) < oo such that (u,a) G T if and only if 
a G (a(u),uj(u)] 

(iv) for any m€T and j G N such that uj G T, a(uj) G (a(u),iv(u)). 

(v) for any u G T and i,j G N such that ui, uj G T, 

M 3 => oi(ui) / a{uj). 
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For any u € T, a(u) is the birth level of u, uj(u) its death level, and we denote by ((u) its 
lifespan C( u ) := oj(u) — a(u). 

Observe that (hi) implies that if T is not reduced to p then has a positive lifespan C(0)- 
We will always assume that a(0) = 0. 

The (possibly infinite) number of individuals alive at chronological level r is denoted by E T 

H r = Card{u € T : ct{y) < r < uj(v)} = Card{x € T : P2(x) = r} < oo, 
and (£ T ; r > 0) is usually called the width process. 








uj{u) 
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Figure 1: A representation of a chronological tree T with finite discrete part T. Horizontal axis 
has no interpretation, but horizontal arrows indicate filiation; vertical axis indicates chrono- 
logical levels. Three elements of T are shown, its root 0, a typical individual v and her mother 
u. Various elements of T are shown, its root p, the death level uj(u) of u, as well as the birth 
levels a(u) and a(v). 

Genealogical and metric structures. A chronological tree can naturally be equipped with 
the following genealogical structure and metric. For any x, y € T such that x = (u, a) and 
y = (v, r), we will say that x is an ancestor of y, and write x -< y as for discrete trees, if u -< v 
and: 

• if u = v, then a < r 

• if u ^ v, then a < a(uj), where j is the unique integer such that uj -< v. 
For y = (v,t), the segment [p,y] is the set of ancestors of y, that is 

[p,y] := {x G T : x -< y} 

= {(v , a) : a{v) < a < r} U {(u, a) : 3k, u = v\k, a(v\k) < a < a{v\k + 1)}. 

For any x, y G T, it is not difficult to see that there is a unique existence point z £ T such that 
[p, x] (~l [p, y] = [p, z]. This point is the point of highest level in T such that z -< x and z -< y. 
In particular, notice that p\(z) = p\(x) /\pi(y) (i.e. p\{z) is the mrca of pi(x) and p\{y))- The 
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level p 2 (z) is called the coalescence level of x and y, and z the coalescence point (or most recent 
common ancestor) of x and y, denoted as for discrete trees by z = x A y. The segment [x, y] is 
then defined as 

[x,y] ■= [p,x] U [p,y]\[p,x Ay[, 

where a reversed bracket means that the corresponding extremity is excluded. A natural 
distance d on T can readily be defined as 

d(x,y) :=p 2 {x) + p 2 {y) -2p 2 (x Ay). 

Note that p 2 is also the distance to the root. 

Degree. The degree of a point x G T, i.e. the number of connected components of T\{x}, 
can be equal to 1, 2 or 3. Apart from the root p, points of degree 1 are called death points or 
leaves and are those x = (u, a) such that a = u(u). Points of degree 2 are called simple points. 
Points of degree 3 are called birth points or branching points and are those x = (u, a) such that 
a = a(uj) for some integer j < K u . For example, if x is not an ancestor of y and y is not an 
ancestor of x, then x A y is a branching point. 

Grafting. Let T, T' be two chronological trees and x = (u, a) G T a point of degree 2. For 
any positive integer i, we denote by g(T', T, x, i) the tree obtained by grafting T 1 on T at x, as 
descending from ui. More precisely, denote by T the tree obtained from T by renaming points 
(ukw,r) as (u(k + 1)w,t), for all k > i and finite integer words w such that (ukw,r) G T. 
Then g(T', T, x, i) is given by 

g(T', T, x, i) := f U {(uiw, a + r) : (u>, r) € T'}. 

Observe that g(T',T,x,i) indeed is a chronological tree. 

Planar embedding. The trees we consider can be regarded as plane trees satisfying the 
rule: 'edges always grow to the right'. For any x G T, we denote by 0(x) the descendance of 
x, that is, the subset of T containing all z G T such that x -< z. The descendance of x can be 
split out into its /(eft)-descendance 6i(x) and r(ight)-descendance 9 r (x). Their definitions are 
as follows: if x is not a branching point, 0i(x) = 9(x) and 9 r (x) = 0; if x = (n, a) is a branching 
point, then a = ot(uj) for some integer j < K u and 

9i(x) := (J (9(u,cr + e) and (9 r (x) := {x} U |J 0(ty',cr + £). 

£>0 £>0 

Actually, since by definition x belongs to its descendance 9(x), it has to belong to either 9i(x) 
or 9 r {x). In the case when x is a branching point, the most convenient convention (which we 
adopt) is that x G 9 r (x). 

Then for any x G T, the complement in T of [p, x] U 9(x) can be partitioned into two forests 
that we call its left-hand component L(x) and its right-hand component R(x) as follows. 

L(x):= (J Biiy) and R(x) := \J 9 r (y). 

y: x&e r {y) y: xeOi(y) 
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For ieT, the branching points that are ancestor points of x belong to one of the following 
sets 

G(x) := {y -A x : x G 6 r (y)} or D(x) := {y A x : x r (y) ^ 0}. 

Elements of G(x) (resp.) are called left (resp. right) branching ancestor points of x, in short 
Ibap (resp. rbap). 

Linear order '<'. Let x, y G T. Assume that x A y g" {x, y}. Then either y G # r (x A y) (and 
then x G 9i(x Ay)) or x G 6 r (x Ay) (and then y G 6i(x Ay)). As a consequence, either x G L(y) 
and y G R(x), or y G L(x) and x G R(y), so that the relation '<' defined on T as follows is a 
total order, or linear order on T (whereas '-<' only defines a partial order): 

x<y<^[y~<xorx£ L(y)] <3> [y A x or y G . 

Note that the genealogical order and the linear order are opposite on a segment [p,x]. This 
will have important consequences when our trees are random. 

It is also important to notice that if T is not reduced to the root, then for any x G T, 

(0,w(0)) < x < (0,q 




Figure 2: Three points x, y, z in a finite chronological tree, satisfying y ~< x and x < y < z. 



Lebesgue measure. The Borel cr-field of a chronological tree T can be defined as the cr-field 
generated by segments. Defining the Lebesgue measure of a segment as the distance between 
its extremities, Caratheodory's theorem ensures the existence of Lebesgue measure, say A, on 
the Borel sets of T. 

We will most of the time call length the measure A(S) of a Borel subset S of T. Note that 
the total length of the tree A(T) is the sum of all lifespans 

A(T) = £C(«)<°°- 
ueT 



8 



We define the truncation of S at level r as the subset of points of § whose existence level is 
lower than r 

C T (S) :={xeS:p 2 (x)<T}. 

We will say that § 

• has locally finite length if for any finite level r, C T (§) has finite length 

• is finite if it has finite discrete part pi(§) 

• is locally finite if for any finite level r, C T (S) is finite. 
Recall that H r is the number of individuals alive at level r. 
Proposition 2.1 For any t > 0, 

' T ~ a da = A(C r (T)). (2) 



o 

7/T /ias locally finite length, then 

• H T < oo /or Lebesgue-a.e. r 

• T /ias finite length iff3 T = for all sufficiently large r. 
Proof. First check that 

A(C T (T)) = ^((^)Ar)-(a(u)Ar)). 

Next write H CT as 

"(J = ^ ] 1{q(«)<o-<w(m)} i 

and use Fubini's theorem to get 

rr rr 

/ E a da = V / dcr l{ a ( u )« 7 < a; ( u )} = V ((^(ti) At) - (a(u) A r)) , 

which ends the proof of ([2]). 

The first item of the list is a mere re-statement of ([2]). To see the direct sense of the 
equivalence stated in the second item, pick any individual u in the discrete part T of T. If 
n = \u\, and u k = u\k for k = 0, 1, . . . , n, then since a{uk) < uj(uk-i), we get by induction 

n 

lo(u) = C(u n ) + a(u n ) < ((u n ) + io{u n ^i) = ((u n ) + C(«n-l) + oi{u n -i) < • < ^ ((u k ). 

k=0 

But when T has finite length E, Y2k=o C( u k) < J2 V <=T^( V ) = ^> so * na ^ u ( u ) < ^- Then all 
individuals are dead at time £, so that H r = for any t > £. 

The converse is elementary. Indeed, if there is To such that H T = for all r > t$, then 
C T (T) = T for all r > to- So if T has locally finite length, it has finite length. □ 

Finally, we indicate how to locally close splitting trees that have locally finite length. 
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Local closure. Let dT denote the (possibly empty) local boundary of the discrete tree T 

dT := {infinite sequences u : u\n 6 T for all n and (a(u\n)) n is bounded}. 

We will sometimes write T:=TU dT for the local closure of T. 

Now assume that T has locally finite length. Then for any u £ dT, one has hin^^oo ((u\n) = 
0, so we can define 

v{u) := lim j a(u\n) = lim uj{u\n). 

n|oo ntoo 

This allows to define the local boundary of T as 

dT : = {(u,u(u)) : u G dT}. 

We point out that this closure is local in the sense that all points in dT are at finite distance 
from the root. Also note that the set dT is not necessarily countable (examples of uncountable 
boundaries will be seen later when trees are random). Points in dT can be thought of as points 
with zero lifespan and infinite height in the discrete genealogy. They have no descendance and 
degree 1, so to distinguish them from other leaves, we will call them limit leaves. Note that 
our terminology is a little bit abusive, since the local boundary of T, taken in the usual sense 
(boundary associated to the metric d), should also comprise actual leaves, as well as the root. 
Also note that limit leaves are their own mother in the discrete genealogy. 

All other properties of T (genealogy, metric, order,...) trivially extend to its local closure 
T := T U dT, and details need not be written down. In particular, Lebesgue measure extends 
to T: its Borel cr-field is the Borel cr-field of T completed with all sets A U B, where A is a 
Borel set of T and B is any subset of dT, with \(A U B) := \(A). Note that A(OT) = 0. 

2.3 Spectrally positive Levy processes 

All results stated in this subsection are well-known and can be found in [3J. 

We denote by (Yt;t > 0) a real-valued Levy process (i.e., a process with independent 
and homogeneous increments, and a.s. cadlag paths) with no negative jumps, and by P x its 
distribution conditional on Yq = x. Its Laplace exponent ip is defined by 

E (exp(-\Y t )) = exp(#(A)), 

and is specified by the Levy-Khinchin formula 

/•oo 

ip(\) = a\ + f3X 2 + / (e" Ar - 1 + Xrl r<1 )A(dr) A > 0, (3) 
J o 

where a G R, f3 > denotes the Gaussian coefficient, and the Levy measure A is a u-finite 
measure on (0, oo] such that J^°(r 2 A l)A(dr) < oo. We will sometimes assume that q := 
A({+oo}) = — ip(0) is positive. This amounts to killing the process at rate q. 

The paths of Y have finite variation a.s. if and only if (3 = and f rA(dr) < oo. Otherwise 
the paths of Y have infinite variation a.s. 

When Y has increasing paths a.s., it is called a subordinator. In that case, ^(A) < for 
any positive A, and we will prefer to define its Laplace exponent as —ip. Since a subordinator 
has finite variation, its Laplace exponent can be written as 

;>oo 

-V>(A) = dA+/ (1 -e- Xr )A(dr) A>0, 
Jo 
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where d > is called the drift coefficient. 

Next assume that Y is not a subordinator. Then, since ip is convex and ultimately positive, 

lim ip(\) = +00. 

A^oo 

Denote by 77 the largest root of ip. If ip(0) < (case when A({+oo}) 7^ 0), then 77 > is the 
unique root of ip. If ip(0) = 0, then either ip'(0 + ) < and 77 > 0, so that ip has exactly two 
roots (0 and 77), or ip'(0 + ) > and ip has a unique root 77 = 0. The inverse of ip on [77, 00) is 
denoted by (p : [0, 00) — > [77, 00) and has in particular 0(0) = 77. 

We write = inf{t > : Y t e A} for the first entrance time of Y in a Borel set A of R, 
and T y for T^}. It is known that 

E (e- qT -*) = e"^ q)x q > 0, x > 0. 

In particular, Po(^-x < 00) = e~ vx . 

Finally, there exists a unique continuous function W : [0, +00) — > [0, +00), with Laplace 
transform 

/ e~ A:r W(x)(ix = — — A > 77, 

such that for any < x < a, 

Px(T < r (O)+00) ) = VF(a - x)/VF(a). (4) 
The function W is strictly increasing and called the scale function. 



2.4 Jirina processes 

We call Jirina process a branching process in discrete time and continuous state-space. Specif- 
ically, a Jirina process is a time-homogeneous Markov chain (Z n ; n > 0) with values in [0, +00] 
satisfying the branching property (w.r.t. initial condition). Writing Z n (x) for the value at 
generation n of the Jirina process starting from Zq = x G [0, +00), the branching property 
implies that for each integer n, (Z n (x) ;x > 0) has i.i.d. nonnegative increments. In particular, 
(Z\(x) ; x > 0) is a subordinator, that we prefer to denote S. Let F be the Laplace exponent 
of S, d its drift coefficient and A its Levy measure. 

By the (homogeneous) Markov property, there are i.i.d. subordinators (S' n ) n >i distributed 
as S, such that, conditional on Zq, Zi, . . . , Z n , 

Z n +i = S n+ i o z n . 

In particular, by Bochner's subordination, the process x 1— > Z n (x) is a subordinator with 
Laplace exponent F n the n-th iterate of F, so that 

E x (exp(-XZ n )) = exp(-xF„(A)) X,x > 0. 

We say that Z is & Jirina process with branching mechanism F. We will sometimes write 
m := F'(0 + ) = d + J °° rA(dr) < 00. If F(0) = 0, the Jirina process is said subcritical, critical 
or supercritical according to whether m < 1, = 1 or > 1. If F(0) > (i.e. A({+cx)}) > 0), it 
is said supercritical. 
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We keep the same symbol rj as that used in the previous subsection, for the following 
quantity 

rj := sup{A > : F(X) = A} < oo, 

because we have in mind that the Levy processes we will consider later on will have Laplace 
exponent A i — s- A — F(X). Last, define the event of extinction, denoted Ext, as 

Ext := { lim Z n = 0}. 



Proposition 2.2 Let Z be a Jirina process with branching mechanism F . Then 

P^Ext) = e - "* x > 0. 

Furthermore, 



{ Z n < oo} = Ext a.s. 



n>0 

Notice that if A is infinite or d > 0, then Z n > a.s. at each generation n, even on extinction. 
On the other hand, if S is a compound Poisson process, then the Borel-Cantelli lemma ensures 
that on the event of extinction, Z n = for all sufficiently large n. 

Proof. Recall that F is concave increasing, so that F{x) > x for x £ [0, rf\ and F(x) < x for 
x > rj. This entails the convergence to rj of the sequence (F n (X)) n for any initial value A > 0. 
As a consequence, 

lim E x (exp(-AZ n )) = e"^ x, X > 0. 
This last convergence implies that P x (Ext) < exp(— rjx), because 

e'^ > lim E I (exp(-AZ n ),Ext) = P x (Ext) x, X > 0. 

n^oo 

Now for any real number A > 0, define by induction the sequence (v n (X)) n as vo(X) = X and 

v n+1 (X) = X + F(v n (X)) n>0. 
Then an elementary induction argument shows that 

E x ^exp ^— A ^2 ^ = exp(-xw n (A)) n > 0. 

On the one hand, by definition, v n+ \ = G\{v n ), where G\ : y i— > A + F(y) is concave increasing 
and has G\(0) > 0. This implies that G\ has at least one fix point, and that its largest fix 
point, say 0(A), satisfying 

4>{X) — F o cf)(X) = X A>0, (5) 

is also the limit of the sequence (t> n (A)). On the other hand, Lebesgue convergence theorem 
gives the Laplace transform of 



T:=^Z,€(0,« 



i=0 
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as 

E x (exp(-\T)) = exp(-x0(A)). (6) 

As a consequence, T < oo with probability exp(— x0(O)). But thanks to ([5]), </>(0) is the largest 
fix point of F, so (f)(0) = r]. Now since {T < 00} C Ext, we get 

e~ vx > Pa, (Ext) > P X (T < 00) = e"^, 

which implies that {T < 00} and Ext have the same probability and hence coincide a.s. □ 

Here the proof is over, but we want to point out that according to ([5]), <p is the inverse of 
A 1— > A — -F(A), which itself is the Laplace exponent of the spectrally positive Levy process 
t 1— > Y t = St — t. Referring to the last subsection on Levy processes, equation ([6]) thus shows 
that T has the same law under P^ as the first hitting time of by the Levy process Y started 
at x. One of the goals of this paper is to shed some light on this relationship. 

3 The exploration process 
3.1 Definition 

In most examples of M-trees, the latter are defined from their contour [H [21 \TE[ I20|. which is a 
real function coding the genealogy, whereas here, we do the opposite (but see also [8]). 

Hereafter, T denotes a chronological tree with finite total length £ := A(T), The real interval 
[0,£] is equipped with its Borel cr-field and Lebesgue measure, which we denote by 'Leb'. 

For any ifT, set 

S(sc) := {y G T : y < x}. 

Since §(x) \ {x} is the union of segments of the form ]z, y] where z = y A x and y ranges over 
the leaves < x, and since the leaves of T are in one-to-one correspondence with T, which is at 
most countable, then §(x) is a Borel subset of T. It is then standard measure theory to prove 
that the real- valued mapping (p : T 1— > [0, £] defined by 

(p(x) := A(S(x)) x G T, 

is measurable and pushes forward Lebesgue measure on T to Lebesgue measure on [0, £]. For 
the sake of conciseness, we will say that ip preserves the (Lebesgue) measure. By construction, 
it also preserves the order. Note that </?(0,u;(0)) = and ip(p) = I. 

Let us show that ip is one-to-one and onto. First, because <p preserves the order, y?(x) = (p(y) 
implies that x < y and y < x, so that x = y, which proves that (p is one-to-one. Second, let 
t G [0,£], and assume that t G D := [0,1] \ ^(T). We want to show that then t is the image 
by (p of a limit leaf. As a first step, notice that because ip preserves the measure, D has zero 
Lebesgue measure. As a second step, we define for each v G T 

g(v) := (p(v, uj(v)) and d[y ) := c^(u, a(v)), 

where u denotes v 's mother. Then observe that for each x G T such that v -< pi(x), one has 

g(v) < ip(x) < d(v), 
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so that in particular, either p\(x) = or there is v with \v\ = 1 such that (p(x) G [g(v),d(v)). 
Assume that there is no v with height 1 such that t G [g(v),d(v)), and let 

U\ := {v G T : \v\ = l,g{v) < d(v) < t} and U 2 := {v G T : \v\ = l,t< g{y) < d(v)}, 

so that U\ and U2 form a partition of the first generation in T. Then we set 

o\ := inf a(v) and a 2 := sup a(v). 

veUi veU , 2 

For each i = 1, 2, if (0, Uj) is a branching point, then there is an integer such that a, = a(k) and 
we put Xi := (k,co(k)), otherwise we put Xi := (0, <7j). Next observe that supj, eWl d(t>) = <p(xi) 
and inf„ G ^ 2 (7(f) = ip(x2), so that 

<^(xi) < t < <p(x 2 ). 

Actually, since by assumption t is not in f(T), we get f(xi) < t < <p(x2), so in particular 
x\ 7^ X2 and o\ > 02- Now for any x G T such that xi < x < X2, and any pair t^) G x 
we have 

g(vi) < d(vi) < <p{x) < g(v 2 ) < d(v 2 ), 

so that pi(x) = 0, and x = (0, a) for some a G (02, 01). Then it is easily seen that <£>(x) = 
9?(xi) + A([x,xi]) = <p(xi) + <7i — <7, and that this equality still holds for x = X2- Therefore 
we get a contradiction, since we can always find a such that t = <^(0, a). In conclusion, there 
must be v with height 1 such that t G (g(v), d(v)) (by assumption, t 7^ g(v))- 

Repeating this argument to the subtree descending from v, and iterating, we get the exis- 
tence of a unique infinite sequence u G dT such that t G (g(u\n), d(u\n)) for all n. These inter- 
vals are nested, and because T has finite length, they have vanishing diameter (d(u\n) — g(u\n) 
is the length of the subtree descending from u\ri), so we get, as announced, t = ip(u, v(u)) for 
some u G dT. 

As a conclusion, 99 is a bijection from T onto [0, £] which preserves the order and the measure. 

Conversely, let ip be any order-preserving and measure-preserving bijection from T onto [0, £]. 
First, since ip is order-preserving, ip~ 1 ([0,tp(x)}) = §(x). Second, since ip is measure-preserving, 
A(S(x)) = A(S(x)) = Leb([0, tp(x)]) = ip(x). This can be recorded in the following statement. 

Theorem 3.1 The mapping ip is the unique order-preserving and measure-preserving bijection 
from T onto the real interval [0, £] . 

Definition 3.2 The process (c ) 5~ 1 (t);t G [0,£]) is called the exploration process. Its second 
projection will be denoted by (X t ; t G [0, £]) and called JCCP, standing for jumping chronological 
contour process. 

Theorem 3.3 The exploration process is cadlag (w.r.t. the distance d on T), and for any 
t G [0,£], t is a jump time iff if" 1 {t) is a leaf (v,oj(v)) ofT. In that case, Lp~ l {t— ) = (u,a(v)), 
where u is v's mother in the discrete genealogy. 

As a consequence, the JCCP (X t ;t G [0,£]) is also cadlag and the size of each of its jumps 
is the lifespan of one individual. In addition, 

X t = -t+ Y, C(«) 0<t<£. (7) 

ip(v,u(v))<t 
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The JCCP is a cadlag function taking the values of all levels of all points in T, once and once 
only, starting at the death level of the ancestor. In the finite case, it follows this rule: when 
the visit of an individual v with lifespan (a(v),uj(v)] begins, the value of the JCCP is uj{v). 
The JCCP then visits lower chronological levels of v's lifespan at constant speed —1. If v has 
no child, then this visit lasts exactly the lifespan £(v) of v; if v has at least one child, then 
the visit is interrupted each time a birth level of one of u's daughters, say w, is encountered 
(youngest child first since the visit started at the death level). At this point, the JCCP jumps 
from a(w) to uj(w) and starts the visit of the existence levels of w. Since the tree has finite 
length, the visit of v has to terminate: it does so at the chronological level a(v) and continues 
the exploration of the existence levels of v's mother, at the level where it had been interrupted. 
This procedure then goes on recursively until level is encountered (0 = a(0) = birth level of 
the root). 

A chronological tree and its associated JCCP are represented on Fig. El which, because of 
its size, was moved to page 03j 

Remark 1 In the case when the tree has finite discrete part, the JCCP has another interpre- 
tation f20[ Fig.l p. 230] in terms of queues. Each jump A t is interpreted as a customer of a 
one-server queue arrived at time t with a load Aj. This server treats the customers' loads at 
constant speed 1 and has priority LIFO (last in - first out). The tree structure is derived from 
the following rule: each customer is the mother of all customers who interrupted her while she 
was being served. A natural ranking of siblings (customers who interrupt the same service) is 
the order of their arrivals. The value Xt of the JCCP is the remaining load in the server at 
time t. 

Now we make two statements that the proof of Theorem 13.31 will require. 

Claim 1. For any x, y S T such that x -< y, we have 

d(x,y) < tp(x) - tp(y). 

Claim 2. Let x G T and (x n )n>o a sequence of points of T converging to x, such that one of 
the following conditions holds 

(a) xo -< x\ < X2 -< • • • -< x 

(b) x -<•••-< X2 -< x\ -< xq and either x is not a branching point, or it is one but x n G 9 r (x) 
for all n. 

Then lim^oo ip(x n ) = ip(x). 

Proof of Claim 1. Recall that d{x,y) = X([x, y\) and notice that E>(x) D S(y) U [x,y], so 
that, taking the measure of each side, 

A(S(x)) > A(S(y)) + A([x, y}) - A(S(y) D [x, y]). 

But S(y) n [x,y] = {y} and \({y}) = 0, so we get d(x,y) = X([x,y}) < A(S(x)) - A(S(j/)) = 
(p(x)-tp(y). □ 
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Proof of Claim 2. Since the proofs in both cases (a) and (b) are very similar, we only write 
it in case (a). Set t := <p(x) and s n := <p(x n ). The genealogical ordering of the sequence induces 
the following ranking xq > x\ > %2 > ■ ■ ■ > x, so that (s n ) is a nonincreasing sequence whose 
limit s := linij^oo s n has s > t. Let us suppose that s > t, and define y := + 1)/2). In 

particular, x < y < x n and y ^ x,x n for all n. Now let z := x Ay. First, because x < y and 
x 7^ y, we have z ^ x. Second, observe that for any n, z and x n are on the segment [p, x] so 
either x n ~< z or z ~< x n ; but since x < y < x n and x n -< x, we deduce that z < x n , so that 
x n -< 2. In conclusion, x n -< z -< x, so that d(x n , x) = d(x n , z) + d(z, x). This brings about the 
contradiction, since d(z,x) > 0, whereas d(x n ,x) vanishes. □ 

Proof of Theorem 13.31 Let x £ OT and i = (p(x). Recall from the beginning of this section 
that there is an infinite sequence u € dT such that for any n, t € (g(u\n), d(u\n)). Also notice 
that the intervals (g(u\n),d(u\n)) form a sequence of nested intervals decreasing to t. But for 
any r, s in this interval, (/9 _1 (r) and ip~ 1 (s) are in the chronological subtree T n descending from 
u\n, so that the distance between those two points in T is at most X(T n ) = d(u\n) — g(u\n), 
which vanishes as n — > oo. This shows that </? -1 is continuous at t, and allows us to discard dT 
in the remainder of the proof. 

Let us show that ip~ 1 is right-continuous. Let t £ [0,£), and write y = ip~ l {t). Next let 
(t n )n>o be a decreasing sequence converging to t, and write y n = (p~ l {t n ). For any z £]p,y[, if 
2/n A y $]z,y[, then y n is explored after z, that is, tp(y n ) — > ^(y)- Since by assumption 
v(yn) converges to v?(y)> we conclude that y n Ay €]z,y[ for all sufficiently large n. This yields 
the convergence of y n A y to y, and since y® Ay ^ y\ Ay < y2 Ay ^ ■ ■ ■ ^ y, we can apply Claim 
2 (a) to the sequence (y n A y) n , which gives 

lim </?(y n Ay) = t. 

n—>oo 

Next, since y n Ay< y n , we can apply Claim 1, which gives d(y n A y, y n ) < ip(y n Ay) - </?(y„). 
But both terms in the r.h.s. of the foregoing inequality converge to t, so that d(y n A y, y n ) 
vanishes. We conclude with the triangular inequality d(y n ,y) < d(y n , y n Ay) + d(y n Ay, y) which 
implies that d(y n ,y) vanishes. In other words, for any decreasing sequence (t n ) converging to 
t, (£> -1 (t n ) converges to i^" 1 ^), that is, ip~ l is right-continuous. 

Next, we prove that ip^ 1 has left-limits that can be characterized as in the theorem. Simi- 
larly as previously, let t € (0,£], and write (v,t) = y = ip~ l {t). This time, let (t n ) n >o be an 
increasing sequence converging to t, and write y n = (/9 _1 (t n ). 

First, assume that y is not a leaf. If y is not a branching point either, then the proof that 
(y n ) converges to y can easily be adapted from that for the right-continuity, but appealing to 
Claim 2 (b) rather than (a). The same argument still applies if y is a branching point, but to 
be allowed to appeal to Claim 2 (b), one first has to prove that y n £ r (y) for all sufficiently 
large n. Let us check that. For any n and any z £ 9 r (y), if y„ #r(y) then y n < z < y, so 
that (p(y n ) < ^(z) < f(y)- But by assumption (f(y n )) converges to (f(y), which proves that 
y n £ 9 r (y) for all sufficiently large n. 

Second, more interestingly, assume that y is a leaf and y £ T, so that r = lu(v). Write 
u for the mother of v, and set z := (u,a(v)), and z' := (u,u>(u)). Since y n < y, y n £" 6 r (z). 
Moreover, applying the same argument as in the end of the last paragraph, we find that for 
all sufficiently large n, y n £ 0i(z), or otherwise said, y n A z' £ [z, z'[. Actually, this same 
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reasoning can be applied to any z" E [z, z'[, namely y n A z' E [z, z"[ for all sufficiently large n, 
which proves that (y n A z') converges to z. Next, let z n be the midpoint of the segment [z, y n \. 
Since y n E &i(z) for n large enough, we have z -< z n -< y n , so that (p(y n ) < ^(^n) an d because 
y E R{z n ), ip(z n ) < t. The immediate consequence is that 

lim (p(z n ) = t. 
n—*oo 

Now applying Claim 1 to z n -< y n , we get d(z n ,y n ) < tp(z n ) — (p(y n ), and since both terms 
in the r.h.s. converge to t, we deduce that d(z n ,y n ) vanishes. But by definition of z n , 
d(z, y n ) = 2d(z n ,y n ), which shows that (y n ) converges to z. In other words (p~ l (t— ) = (u, a(v)). 

The proof now focuses on the JCCP X = (Xt\t E [0, £[). First, X = p 2 o tp^ 1 is cadlag 
as a mere consequence of the fact that ip -1 is cadlag and p 2 is continuous (whereas p\ is not, 
though). Now for any t E (0, 1), t is a jump time of X only if it is a jump time of y? -1 . When 
it is so, we know that <^ -1 (t) = (v,lu(v)) for some individual v, and that c/9 — 1 (t — ) = (u,a(v)), 
where u is v's mother. This shows that Xt = ^(v) and Xt- = a(v), so that the jump size at 
time t, AX t = X t — X t _, is equal to the lifespan C(-u) of v. 

Now we prove ([7]) . As said in the beginning of the proof, we can assume that x = ip~ l (t) E T. 
We will need the following notation: for each point y = (v,t) E T, we will set y := (u,a(v)), 
where u is u's mother. Notice that since y ~< y, we always have y < y. Then S(x) can be 
written as the union of C(x) and U(x), where 

C(x) := {y E T : y < y < x] and U(x) := {y E T : y < x < y}. 

The only intersection between these subsets is x, and only if x is a branching point. We will use 
twice the following observation: if (v,t) E S(x), then for any a E [r, u>(v)], since (v,t) -< (v,a), 
we have (v,a) < (v,t) < x, so that (v,a) E §(x). The first consequence of this observation is 
that for any y E C(x), with v = pi(y), (v,a) E C(x) for any a E (a(v ), u>(v)]. We write C(x) 
the set of such vertices of the discrete tree. The second consequence is that for any y E U(x), 
with y ^ x and v = p\{y), there is r E (a(v),u>(v)) such that 

(v,a) E U(x) -v4> <7 E (t,u(v)\, 

which implies that (v,a) < x < (v,t). Now thanks to Claim 1, unless (v,t) = x, z = (v,t) 
is a branching point, (v,a) E 0i(z) and x E r (z). In particular, if w denotes p\{x), we have 
v -< w in the discrete tree, so that pi(U(x)) is the set of ancestors of w, which has cardinality 
n + 1, where n = |to|, and, thanks to the last display, 

n-1 

U(x) = (w x [p2(^)i «(«)]) U ((u)|A;) x (a(w\k + l),u(w\k)]) . 

k=0 

Taking the Lebesgue measure, we get 

n-1 

A(U(x)) = u(w) -p 2 (x) + ^2(u(w\k) - a(w\k + 1)) 

k=0 

On the other hand, 

A(C(x)) = ^ C(«), 

ueC(a:) 
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so that, recalling that the intersection of C(x) and U(x) has zero Lebesgue measure, and that 
their union equals we deduce 

A(S(x))=[ ((v)\ - (p 2 (x)-a(w)+^(a(w\k + l)-a(w\k))Y 

\v:(v,w(v))<x J V fc=0 / 

Now notice that all terms in the second expression in the r.h.s. cancel out to p 2 (x), which 
yields 

p 2 (x) = -tp(x) + ^2 CO)- 

v.(v,w(v))<x 

Writing t = tp(x) yields ©. □ 
3.2 Properties of the JCCP 

Actually, the chronological tree itself can be recovered from its JCCP (modulo labelling of sib- 
lings). In the next two statements, we provide some useful applications of this correspondence. 
For each t G [0, £}, set 

i : = sup{s < t : X s < X t ] V < t < I. 



Theorem 3.4 Let x = (u, a) and y = (v, r) denote any two points in T, and set s = (f(x) and 
t = (p(y). Then the following hold : 

(i) The first visit to v is t 

(p(v, uj(v)) = t. 

In particular, if y G dT then t = t. If t is a jump time of X, then t = t as well, and the first 
visit to the mother u of v in T is given by 

ip(u,uj(u)) = sup{s < t : X s < X t -}. 

(ii) Ancestry between x and y : 

y ^ x i < s <t 
(Hi) Coalescence level between x and y (assume e.g. s <t) : 

p 2 (x Ay) = inf X r . 

s<r<t 

For any t G [0,£], we define the process X® on [0, t] as 

X® := X t - - X {t _ r) _ re[0,t], 
with the convention that Xq_ = 0. We also set 

H t := Ipiap- 1 ^ 
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the generation, or genealogical height in T, of the individual v = p\ o ip~ 1 (t) visited at time t 
(recall that \v\ denotes the length of the integer word v). 

The following corollary states that if v has finite height Ht, then the record times of are 
exactly those times when each of v's ancestors is visited for the first time by the exploration 
process (which actually holds also when v € dT). It also characterizes the height process 
(Ht; t > 0) of genealogical heights, or height process, as a functional of the path of the JCCP. 

Corollary 3.5 Let y = <p (t), v = p\(y), and the first visit to vj. = v\k (ancestor of v 
belonging to generation k), that is, 

(i) if n := \v\ < oo, then y dT, and one can define recursively the record times of X® by 
s\ = t — t and 

s k+ i = mf{s > s k : X« > X«} k > 0. 

Then 

tk = t- s n _k+i < k < n. 

(ii) in the general case, recall Ht < oo is the genealogical height ofv^T. Then Ht is given 

by 

H t = Card{0 < a < t : X® = sup X®} 

0<r<s 

= Card{0 < s < t : X s _ < inf X r }. (8) 

s<r<t 



Quantities defined in the previous two statements are represented in Fig. [3j 




h t 2 h 1 v 1 

Figure 3: The JCCP of some finite chronological tree with jumps in solid line. Set u = p\OLp~ l (t) 
(resp. v = p\ o ip^it')) the individual visited at time t (resp. t 1 ). The first time when v is 
visited by the exploration process is t'. The first visits to the Ht = 3 ancestors of u are also 
shown. 
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Proof of Theorem 13.41 (i) We first show that when y £ dT, then t = t. Indeed, recall that 
v is an infinite sequence, and that t = lim^oo f t n , where t n = (f(v\n, u>(v\n)). Also thanks to 
Theorem 13. 3| X tn _ = a(v \n), so X tn _ increases to X t - = X t . This shows that i = t. 

As a by-product we also get the displayed equation in the case when y G dT, since then 
there is only one pointwise visit t of v, so the first visit is t and we know that t = t. From now 
on, we can discard the points of dT. 

Set z = (v,u>(v)), where v = p\ o (p~ l (t) is the individual visited at time t (and recall 
y = <p (t)). Observe that z < y because y -< z, so that y?(z) < t, and recall from Theorem 
13.31 that ip{z) is a jump time, with X^ z )_ = a(v) < P2(y) = Xt. Now for any s £ [(p(z),t], set 
x = and a = x A z. Since both a and y belong to [p, z], either a ~< y or y -< a. We will 

suppose that a -< y and find a contradiction. Recall z < x. Either x -< z, so a = x and y < x, 
or x £ R(z), so that x £ R(y), and similarly y < x. But unless y = x, this contradicts s < t. 
As a consequence, y -< a, so that y ^ x, and A t = p 2 (y) < P2(^) = ^s- Summing everything 
up, Xy( z \ < X t and for any s € X s > X t . From these, we deduce that t = ip(z). 

For the second statement, let t be a jump time, and set z = (u, uj(u)), where u is the mother 
of v in the discrete tree T (it exists since ^^{t) is a leaf and so is not in dT). Then from 
Theorem 13.31 we know that y = (v,uj(v)), that is, t = ip(v,uj{v)) = t, that X t - = a(v), and 
similarly X v r z \ = a(u). Because y € R{z), we get ip(z) < t, and because v is a daughter of u, 
we get a(u) < a(v), that is, X v ^_ < X t -. Now for s G](p(z),t[, set x = ip~ 1 (s) and a = x Az. 
We want to show that X s > X t -, which also reads P2(x) > a(v). Since s > ip(z), either x -< z 
or x £ R(z). If x ~< z, then P2(x) > a(v) otherwise we would get x ~< y and then t < s. If 
x £ R(z), then a is a rbap of z with p2(a) > a(v) otherwise we would get x £ R(y) and then 
s > t. But since a -< x and a/i, P2(%) > Piip), and we deduce again that P2(x) > a(v). So 
we get X s > X t - for all s £]ip(z),t[, and the proof ends recalling that X<pi z \ < Xt—. 

(ii) Again set z = (v,u(v)), with v = pi(y) and y = (p~ l (t). Recall that in the first para- 
graph above, we have shown that for any s £ [ip{z),t], y -< x(= 99 _1 (s)). Now thanks to (i), 
t = tp{z), so that for any s E [t, t], y -< x. Conversely, assume that y < x. First notice that 
x < y, so that s < t. Second, we want to show that z < x, so that t < s. Set a = x A z. Then 
either x = a, so that x -< z (and z < x), or a is a branching ancestor point of z. Now since 
y -< x and y -< z, we get y < a, and a G [y, z]. Finally, because pi(y) = p\(z)(= v), there is no 
lbap in the segment [y, z\, so that a is a rbap of z and x £ R(z) (then z < x again). 

(iii) Set / S)t = inf s < r < t X r . If y -< x, then p 2 (^ A y) = p 2 (y) = ^t- Now, by (ii), i < s <t, 
and by (i), ij ( = A t , so that I Sj t = Xj. Thus, we have proved (iii) when x A y = y. Next 
set z = x A y and assume without loss of generality that y G -R(x), so z is a branching point 
(u,a{y)) = (/? _1 (r). Since z < y, t < r . Also we have z -< x, so f < s, which means that 
Xh > X r for any h £ [s,r] C [s,t\. In particular, I s>t > A r = P2(x Ay). For the converse 
inequality, notice that x < a < y, where a = (u,u;(u)) because a G i?(x) and y G i?(a). Also 
by definition, p\{x) Ap\{y) = u and since y G R(x), v is an ancestor individual of pi(y) but 
not of pi(x) (otherwise pi(x) Ap\{y) = v). This ensures that a ^ x, and since x < a < y, 
s < tp(a) < t. But a is a leaf, so thanks to Theorem 13.31 (^^((^(a) - ) = (u,a(v)), and 
Xtp(a)- = a ( v ) = P2{x Ay). Because s < cp(a) < t, we get I s j <P2{x Ay). □ 
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Proof of Corollary 13.51 We show (i) by descending induction on k £ {0, ...,n}. Since 
\v\ = n, v n = v, and by the first statement in Theorem 13.41 (i), t n = tp(v, u>(v)) = t = t — s%, so 
the result holds for k = n. Next let k < n — 1, and assume that t k+ \ = t — s n _ k . Thanks to 
the second statement in Theorem 13.41 (i). since v k is the mother of v k+ \, 

t k = <p(v k ,u(v k )) = sup{s < tk+i : X s < X tk+1 -}. 

Elementary calculations yield 

sup{s < t k+1 : X s < X tfc+1 _} = t - inf{s > t - t k+1 : X® > X® th+ J. 

But since t — t k +\ = s n - k by assumption, we deduce that 

t k = t - mf{s > s n _ fc : X® > X®_ h } = t - s n . k+l . 

As for (ii), if v dT, just use (i) to see that the number of records of J« is the number 
of ancestors of </? — 1 (£) in T, which is also its height. Then check that the records of the past 
supremum of XW are those of the future infimum jW of X on [0,t], that is, 

= inf X r 0<s<t. 

s<r<t 

Now if v G <9T, then on the one hand v is infinite so by definition Ht = oo; on the other hand, 
apply (ii) to v\n + 1: write t n = ip(v\n,a(v\n + 1)) and recall that X s > X tn for s G [t^ , t] , so 
that the future infimum Jw of X on [0, i] has at least n records, and let n — » oo. □ 

4 Splitting trees 

In this section, we consider random chronological trees, called splitting trees, and whose width 
process, in the locally finite case, is a binary, homogeneous Crump-Mode-Jagers process. 

4.1 Definition 

A splitting tree is a random chronological tree characterized by a u-finite measure A on (0, oo] 
called the lifespan measure, satisfying 

/ (r A l)A(dr) < oo. 

■/(0,oo] 

Let P x denote the law of a splitting tree starting with one ancestor individual having deter- 
ministic lifetime (0, x], where x is only allowed to equal oo when A({+oo}) > 0. 

We give a recursive characterization of the family of probability measures P = (P x ) x >o 
as follows. Recall from the Preliminaries on chronological trees that g(T',T, x,i) is the tree 
obtained by grafting T' on T at descending from p\(x)i. 

Let {on, Ci)i>i be the atoms of a Poisson measure on (0, x) x (0, +oo] with intensity measure 
Leb(g> A (where 'Leb' stands for Lebesgue measure). Then P is the unique family of probability 
measures on chronological trees T satisfying 

T= (J <?(T n ,0x(O,x),(0, 

oi n ),n), 

n>l 
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where, conditionally on the Poisson measure, the (T n ) are independent splitting trees, and for 
each integer n, conditional on £ n = £, T n has law P^. 

In other words, for each individual v of the tree, conditional on a(v) and lo(v), the pairs 
(a(vi), C(vi))i>\ made of the birth levels and lifespans of v's offspring are the atoms of a Poisson 
measure on (a(y), u)(v)) x (0,+oo] with intensity measure Leb <8> A. In addition, conditionally 
on this Poisson measure, descending subtrees issued from these offspring are independent. 

Here, we have to say a word about the order of siblings. First, note that the right-hand 
side of the last display can also be understood as the sequential grafting of trees on x (0, x) 
as n increases. Second, as a referee pointed out, we need to specify how we label the atoms 
(«i ; Ci)i>i to characterize P. In the case when A({+co}) = 0, all lifespans are finite a.s., and 
we can assume that atoms can be ranked in such a way that for any two distinct integers i, j, 

i < j <=4> Ci > Cj or (Ci = Cj and cq < a,). (9) 

Note that this order carries over to the case when A({+oo}) > but x < °°- I n the case when 
A({+oo}) > and x = °°j we fi x a bijection h from N onto N x Z + , and we proceed into three 
steps. First, all atoms (a, C) with infinite lifespan (second marginal) are ranked in increasing 
order of their birth date (first marginal). Then we call i-th cluster of atoms, the subset of 
those atoms with finite lifespan whose birth date is between the birth dates of the (i — l)-th 
and i-th atoms with infinite lifespan (birth date of 0-th atom is 0). Second, for each cluster 
separately, we can (and do) rank the atoms of the cluster in the order defined in ([9]). Third 
and last, writing h(n) = (hi(n),h.2(n)), we label atoms altogether in such a way that 

(on>Cn) has rank hvin) in the hi(n)-th cluster, (10) 

where the atom with rank in the i-th cluster is the i-th atom with infinite lifespan itself. 

Check that when A is finite, the splitting tree is locally finite a.s. The terminology 'split- 
ting trees 1 is usually restricted to this case |llj . We will say a few more words about it in the 
next subsection. 

Two branching processes. Recall that H r denotes the number of individuals alive at t 
(width of T at level t) : 

E T = Card{t> G T : a{v) < r < u>(v)} < oo. 

Under P x , the width process H is a branching process: allowing several ancestor individuals 
with lifespans equal to x an d i.i.d. descendances, the branching property holds as a function 
of the number of ancestors. Unless A is exponential or a Dirac mass at {oo}, this branching 
process is not Markovian. 

Next set Z n the sum of all lifespans of individuals belonging to generation n : 

Z n := C(v) = \({x£T:\ Pl (x)\=n}). 

v£T:\v\=n 

Under P x , Z is another branching process, which itself is Markovian, as seen in the next 
statement. 
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Theorem 4.1 Under ¥ x , (Z n ;n > 0) is 
mechanism F given by 

F(X) := / (1 
J 

In addition, T has locally finite length a.s. 
and the following events coincide a.s. 

(i) T has finite length 

(ii) lim n ^oo Z n = 

(Hi) E T = for all sufficiently large r. 

Remark 2 The last statement allows to use the same terminology for splitting trees as for 
Jirina processes. For example, recall from the Preliminaries that the event defined a.s. equiv- 
alently by (i), (ii) and (Hi) is called extinction and denoted Ext. So when F(0) = 0, we set 
m := F'(0 + ) = Jq°° rA(dr), and we say that T is subcritical, critical or supercritical according 
whether m < 1, = 1 or > 1. When F(0) ^ (that is, A({+oo}) ^ 0), T is said supercritical. 



a Jirina process starting from x> with branching 

-e- Xr )A(dr) A>0. 
under P, so H r is a.s. finite for Lebesgue-a.e. t, 



Proof. By construction, 

where (Aj; s > 0) is a Poisson point process with intensity Leb <g> A. Since Jq rA(dr) < oo by 
assumption, Z\ has the value at time x OI a subordinator S with Laplace exponent F. Now we 
reason by induction on the generation number n. By definition, conditionally on the knowledge 
of generation n, 

^ +1 = E E Ai s \ 

v:\v\=n s<£(v) 

where the point processes ( ; s > 0) are i.i.d. Poisson point processes with common intensity 
Leb A. It is then a standard property of Poisson measures that, since Z n — X^w|u|=n C(^0j 
can write 



E 

s<Z n 



Zn+1 — / . A", 



where (A"; s > 0) is a Poisson point process with intensity Leb <8> A. Then by induction on n, 
Z n has the law of the n-th composition (in Bochner's sense) of i.i.d. subordinators distributed 
as S, and evaluated at x- This is precisely saying that Z is a Jirina process starting from x 
with branching mechanism F (see Preliminaries). 

To show that T has locally finite length a.s. we are going to use Proposition 12.11 Set 
the number of individuals from generation n living at time r 

E< n ) := Card{x G T : \pi(x)\ = n,p 2 (x) = r}. 
Observe that = 1{ X > T } < oo, and let n > 0. By definition, conditional on (H^; a < r), 



:(n+l) = V" 1 



{ai<r<a 4 +Ci}) 
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where (aj, Ci)i>i are the atoms of a Poisson measure on (0, r) x (0, oo] with intensity measure 
Ho -da A. This entails that the conditional distribution of H^ n+1 ' ) is Poisson with parameter 

[ T daE^ f A(dr). 

JO J[t— <t,oo] 

With the notations A(r) := A([r, oo]) and f n (j) := E x (Hf^), we thus get 

/n+i(r) = [ T daf n (a)A(T-a). (11) 
J o 

Next denote by G n the Laplace transform of f n 

poo 

G n (ji) := / Utye-^dt /i>0, 
Jo 

and assume that the series of generic term G n {ij) converges for some \x. Then Fubini's theorem 
ensures that t h-> f(t) := fn(t) is integrable against i i— > e~ M *. In particular, for any r > 0, 
J T /(i) (it < e^ T J r /(t) (it < e^ T j °° f(t) dt < oo. Now applying Fubini's theorem to 

Ejj = J2 n ■ > we get E X (H CT ) = /(c). A third and last application of Fubini's theorem, along 
with equation ([SJ in Proposition 12.1} yields 

f'T f'T 

E X A(C T (T)) = E x / Eq- da — I f(a) da < oo for any r > 0. 
Jo Jo 

In conclusion, it is sufficient that Y2 n G n (p) < oo for some n to ensure that T has locally finite 
length a.s. (and, actually, has even integrable local length). Now we prove that it is indeed 
the case that Y2 n G n {n) < oo for some \i. Thanks to (jlip . we have 

/'OO i*t roo 

G n+1 (fx)= dte-i* dsf n (s)A(t-s) = G n (fx) dte~^A{t). 
Jo Jo Jo 

Now a straightforward calculation shows that f °° dt e~^ 1 A(t) = F(fi)/fi, so that 

° M= {PIT)" GM M> °' 

and Gq{(jl) = f x e~^ lt dt < oo for any fj, > 0. Now recall from the Preliminaries that there 
is some finite nonnegative real number r\ such that /j, > rj F(p>) < /i, which entails the 
convergence of ^ n G n (/i) for all \i > rj. 

We finish the proof with the a.s. equality of the events defined in (i), (ii) and (iii). Be- 
cause A(T) = Y2 n >o and Z is a Jirina process, we can apply Proposition 12.21 to get the a.s. 
equality of the events defined in (i) and (ii). The a.s. equality between those defined in (i) and 
(iii) comes from an appeal to Proposition 12.11 recalling that T has a.s. locally finite length. □ 



24 



4.2 The finite case 



The Crump-Mode- J agers process, or CMJ process, is sometimes called the general branching 
process (see [26] for a complete overview), because it counts the size of a branching population 
defined under the very general assumption that 'reproduction schemes' of all individuals are 
i.i.d.. Specifically, the reproduction scheme of an individual with given birth date a is the joint 
knowledge of 

• her lifespan £ £ (0, oo] 

• the successive ages (0 <) u\ < 02 <•••(< C) a t which she gives birth 

• the (integer) size £j of the clutch she begot at the i-th birth time a + Uj, for all i > 1. 

It is clear that a branching population with i.i.d. reproduction schemes can be constructed 
recursively starting from a finite number of individuals with given birth dates. 
To stick to the present framework, we further assume that 

1. all clutch sizes (the £j's) are a.s. equal to 1 (binary splitting) 

2. conditional on the lifespan £, the point process (<7j) is a Poisson point process on (0, Q 
with intensity b (homogeneous reproduction scheme) 

3. the common distribution of lifespans is A(-) /b, where A is some positive measure on (0, 00] 
with mass b called the lifespan measure. 

In other words, each individual gives birth at rate b during her lifetime (a, lj], to independent 
copies of herself whose lifespan common distribution is A(-)/b. Check that this definition of a 
splitting tree is exactly the same as that given previously, in the special case when A is finite. 
In that case, the width process (H r ;r > 0), is a homogeneous binary Crump-Mode- J agers 
process. 

Trees satisfying the foregoing assumptions could also be called general binary trees with 
constant birth rate, or homogeneous binary Crump-Mode-Jagers trees, but we will stick to the 
name of ''splitting trees 1 . On the other hand, this terminology is unfortunate because it evokes 
renewing binary fission (Yule tree). 

A third branching process. In the finite case, one can define Z n the number of individuals 
belonging to generation n 



From the definition, for any individual v, the total offspring number of v, conditional on 
C(v) = z, is a Poisson random variable with parameter bz. Then it is easy to see that (Z n \n > 0) 
is a Bienayme-Galton- Watson process started at 1, with offspring generating function / 



and the per capita mean number of offspring can easily be computed to be equal to m (recall 



2-n '■ — 



Card{v € T : \v\ = n}. 





m = 



J °°rA(dr)). 
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Remark for modeling purpose. This birth-death scheme can be seen alternatively as a 
constant birth intensity measure bLeb combined with an age-dependent death intensity measure 
/i given by 

P(C>*) = exp[-M((0 } *))] z>0, 

which forces the equality 

exp[-/i((0,z))] = A(z)/6 z > 0, 
where A(z) = A([z, oo]), and yields the following equation for [i 

„ ((fa) = , > o. 

It has to be emphasized that in survival analysis (used in systems reliability, medical research, 
actuarial science, conservation biology,...), the (density of the) measure \x is called the hazard 
function (probability of failing during time interval dz conditional on survival up to time z), 
whereas the (density of the) measure A/6 is called the failure rate. 

If for some reason one has to proceed the other way round (the birth rate b and some death 
intensity measure n are given), notice that the lifespan measure A is then 

A{dz) = b/j,(dz) exp[-/i((0, z))] z > 0. 

The requirement that fj, has to fulfill for lifespans to be a.s. finite is fi((0, oo)) = oo (correspond- 
ing to A({+oo}) = 0). On the contrary, if \x is the null measure on (0, oo) (i.e. A({+oo}) = 6), 
then the CM J process is called a Yule process (pure-birth process with constant rate b). 



4.3 Law of the JCCP of a splitting tree 

From now on, we consider a splitting tree T with lifespan measure A, and whose ancestor has 
lifespan £(0) = \- Recall that its law is denoted P x . 

From Theorem 14.11 we know that T has locally finite length a.s., and finite length iff 
extinction occurs. In particular, T has a JCCP on Ext, denoted by (X t , t G [0, A(T)]), and any of 
its finite truncations C T (T) has a JCCP as well, denoted by (X t (r) , t G [0, A(C T (T))]). Although 
similar with the notation used in Corollary 13.51 this one has a totally different meaning. 

Also set 

E T := Card{x G T : p 2 (x) = r} r > 0. 

Observe that for any r' > t, the set {x G T : p 2 {x) = r} is equal to {x G C T /(T) : p 2 {x) = r}, 
so that 

E T = Card{£ : X { p = r} < r < r', 
which does not depend on r'. 

Lemma 4.2 For any fixed x> r ; ^x(" r < °°) = 

An important consequence of this lemma is that for any given \ an d r — T \ both X (on 
extinction) and X^ T '^ hit r a finite number of times a.s. under P x . 
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Proof. Assume that there is X) r > such that b X)T := F X (E T = oo) > 0. Then observe that, by 
the definition of splitting trees, % l— ¥ b x>T ls nondecreasing on [0, r] and constant on [r, +00). 
This ensures that b TjT > 0. But using the same kind of argument, as well as translation 
invariance, we get 6 CTi(J > 6 T)T for any a > r and then b x > j(T > b TjT for any ^ > o~ > t. As a 
consequence, we get that for any %' > r, 

roc 

E X' J Q %=oo}^>(x'-*,r>0. 

Now observe that for any Borel set A G T, Leb(p2( j 4)) < M^)- Moreover, H CT = 00 only if 
B„ = 00 or a G £>2(<9T). But on the one hand, thanks to Theorem 14.11 {a : = 00} has 
a.s. zero Lebesgue measure, and on the other hand, <9T, hence p2(dT), also have zero Lebesgue 
measure. This is in contradiction with the last display. □ 

We denote by Y the spectrally positive Levy process t \—> Y t := — t + J2 s <t where (A t , t > 0) 
is a Poisson point process with intensity measure Leb (g) A. In particular, Y is a Levy process 
with finite variation, whose Laplace exponent (see Preliminaries) will be denoted by ip 

POO 

■0(A) := A - F(X) = A - / (1 - exp(-Ar)) A(dr) A > 0. 

The following statement is the fundamental result of this section. It is a little bit surprising at 
first sight, in the sense that, eventhough (^j (t); t > 0) is not Markovian, its second projection 
is. Recall that Ta is the first hitting time of A. 

Theorem 4.3 The law of is characterized by (i); conditional on Ext, the law of X is 
characterized by (ii). 

(i) Define recursively to = 0, and = inf{t > t\ : X^ G {0, r}}. Then under F x , 

the killed paths ej := (X^ t ,0 < t < tj+i — ti), i > 0, form a sequence of i.i.d. excursions, 
distributed as the Levy process Y killed at Tq A 2} Tj+0O ), ending at the first excursion hitting 
before (r, +00). These excursions all start at r, but the first one, which starts at min(x>T). In 
other words, X^ has the law ofY reflected below r and killed upon hitting 0. 

(ii) Under P x (- | Ext), X has the law of the Levy process Y, started at x> conditioned 
on, and killed upon, hitting 0. 

Proof, (i) In what follows, we will stick to the notation ip~ 1 for the exploration process of 
C T (T) and we_define i := A(C T (T)). Then let t G [0,£) and x = <f~ l {t) G C T (T) = C T (T). Set 
u = p\{x) € T and n = \u\ < 00. Set also = u\k. The segment [p,x[ of ancestor points 
of x (except itself) is the union of segments x [a(uk), a(«fc+i)), for < k < n, with the 
convention a(u n+ i) = P2{x) = xj 7 ^ when n < 00. Next denote by 77^ G (a(uk), u(uf.)) the 
birth level of u^i, and by 1^ the set of integers i such that Tf-i < a(ui.+i), that is, u^i is a child 
of Uk born before u^+x- Finally, we denote by T^j the subtree that can be seen as grafted at 
Tk t i, that is, 

Tfc,i := {(v, a) : (u k iv, a + r^) G 9 r ((uk, < k <n,i <E I k . 
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Now, we can write the subtree lp2, t _ t of points of C r (T) still unvisited at time t (which is the 
union of [p,x[ and R(x) in C T (T)) as 



We point out that because truncation preserves the order, the post-t exploration process of 
C T (T) is exactly the exploration process of T p og t _ t . 

Now we work under F x . We call heightwise label transposition any mapping 8 : U — > U, 
for which there are integers k,i,j such that for any y G U, written as y = (v, a) with u = 



where (resp. v**) is obtained from v by setting equal to i (resp. to j), and leaving other 
labels unchanged. Observe that heightwise label transpositions map chronological trees into 
chronological trees. 

For any r > 0, let T r denote the cr-fleld generated by {F(ip~ 1 (s)); s < r,F G 7i\, where 
7i denotes the set of measurable (real) functions on U that are invariant under the action of 
all heightwise label transpositions. This filtration is chosen to contain all past events that 
do not depend upon the choice of sibling labelling. Notice that since P2 £ H, (X s ;s < r) is 
J>-measurable, so it is sufficient to prove that is Markovian w.r.t. the filtration (T r ) r >o- 

Now for all < k < n, let C,k,i be the lifespan of the ancestor of T^i, and let [i^ be the 
random point measure on (a(uk), a{uk+i)) x (0, oo] with atoms ((Tk,i,(k,i),i £ Ik)- 

Claim 1. For any k = 0, ...n, conditional on {a(uk), a(uk+i)}, is a Poisson point mea- 
sure on (a(iifc), a(uk+i)) x (0, oo] with intensity measure Leb <8> A, independent of F t . 

To see this, notice that because i G 1^ ^ r^^ < a(uk+i), the point measure Hk is the restric- 
tion to (a(ufc), a{uk+i)) x (0, oo] of the Poisson point measure, say v^, on (a(uf~), w(«fc)) x (0, oo] 
with intensity measure Leb<8> A. Now we are reasoning conditional on {a(uk), ot(iik+i)}- First, 
the point measure [i^ (conditionally) is a Poisson point measure on (a(v,k) , ce(uk+i)) x (0, oo] 
with intensity measure Leb£g>A. Second, as a consequence of the recursive construction of T, ^ 
depends on Tt only through the knowledge of ((r^j, (k,i),i & Ik)- But since Tt is independent of 
sibling labelling, [i^ depends on Tt only through the knowledge of the atoms of the restriction 
of ffc to [a(uk+i),Lv(uk)) x (0, oo]. As a consequence, the point measure fik is (conditionally) 
independent of Tt- 

Claim 2. Conditional on (a(uj,);0 < k < n), the point measures (/Ufc;0 < k < n) are 
independent. 

This last claim is a mere consequence of the recursive construction of T, and of the fact 
that all points ((ufc,cr);0 < k < n, a G (a(uk), a(v,k+i))) have disjoint descendances. 



T 




{vi, . . -,Vl), 




(v,a) if / < k or v k {i,j 
(v*, a) if I > k and v k = j 
(u**, a) if I > k and v k = i, 



} 
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Now rank the atoms ((r^j, Ofc,i)> < k < n,i £ If.) in the order specified by Q or (fTUj) . 
according to whether A({+oo}) is zero or not. We denote by (rj,Q;j > 1) this relabelled set 
of atoms, as well as Tj the tree corresponding to (Tj,(j) (read j as j(k,i)). For precisely the 
same reason as for Claim 2, we have 

Claim 3. Conditional on > 1), the trees Tj are independent and independent of Tt, 
and for all j > 1, conditional on = z, Tj has law P^. 

Next define 

T':=U 5 (T J ,0x[O,X, (r) [,(0,r i ),i). 
i 

Thanks to Claims 1 and 2, conditional on -Xj 1 ^ = a, and conditional on (a(u&); < A; < n), the 
point measure fj, := X^o<fc<n Mfe * s the superposition of independent Poisson point measures, 
all independent of Tt, all with the same translation invariant intensity Leb <8> A, on disjoint 
sets whose union is equal to (0, a) x (0, oo], up to the discrete set (a(uk);0 < k < n). As a 
consequence, fi is a Poisson point measure with intensity Leb® A on (0, a) x (0, oo], which does 
not depend on {a(ufc +1 ) : < k < n}, and is hence independent of T% conditional on = a. 
Now thanks to Claim 3, conditional on xj: = a, T' is independent of Tt and has law P CT . 
Now observe that any x £ Tpog t _ t has pi(x) = u^iv for some v , and that the mapping from 

TTpost-t to C T (T') that maps u^iv into jv (j = j(k,i)), and leaves P2{x) unchanged is a bijection 
which preserves the linear order and the second projection. As a consequence, the JCCPs of 
both trees are equal. This shows that conditional on x[ T ^ = a, (X^ s ; s > 0) is independent of 
T and has the distribution of (X^; s > 0) under P CT . 

Thus, we have shown that X^ is Markovian, and since t = inf{s > : X^ = 0}, we get 
that (jW;0< s < i) is a Markov process killed upon hitting 0. 

We want to prove that X^ is a Feller process. Let / be a continuous, hence bounded, function 
on [0,t]. 

As a first step, because X^ T ' is cadlag, the dominated convergence theorem ensures the 
right-continuity of t \— ► E x ff^X^jj, which is the first required condition for X^ to be 
Feller. 

As a second step, let x' < X-, an d recall that (0, %') is P x -a.s. a simple point, that is, a 
point with degree 2 in the chronological tree, at which the exploration process is continuous. 
Also notice that <^((0, x')) is the first hitting time T x > of %' by X^ T >. As a consequence, 
the contour process of the (truncation of the) chronological tree issued from x [0, x'] (resp. 
x [x',x]) is {Xrl+vt > °) ( res P- P^ (t) ;0 < t < T x >)), and in addition both contour 

processes are independent, and the former has the law of X^ under ¥ x >. In particular, X^ 
satisfies the Markov property at T x > under P x . Using these observations, we now show that 

X ► E x ^/(A 7 ^)^ is continuous on [0,r], which is the other required condition for X^ to 

be Feller. First, since 2y is the length of the subtree issued from x [x',x] an d truncated at 
t, it is stochastically dominated by A(T) under P x _ x /. Putting together the first statement of 
Theorem 14.11 with the arguments developed in the proof of Proposition 12.21 we see that under 
F x , A(T) = Y^n (='■ T in the proof of Proposition E2]) converges to in probability as x — > 0. 
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This ensures that F x (T x i > e) vanishes as x ~ x' ~* 0. Second, we have the following equality 

E x (f(X^))-E x ,(f(X^)) = J\ x (f(X^)-f { xiZ),T x ,eds) 

+ E x (/(X f (T) )-/(<T>^> e ). 

Recalling that / is bounded, using dominated convergence along with the right-continuity of 
X^ for the first term, we get the continuity of x *—* (/(-^t^)) - ^ ne conclusion (see e.g. 
[25J) is that X^ is a Feller process, and as such, it satisfies the strong Markov property. 

Now fix t' > t > and define A T,T ' as 

and let a T ' T ' be its right inverse 

a/ = inf{s : A T / > t} t>0. 

Then it is clear from our construction that A"( T ') oa T,T ' is the JCCP of C T oC T '(T). But because 
t' > t, C T o C T i = C T , so that 

X (r) = X (t>) oa r,r'_ 

Then define recursively tq = 0, S{ = inf{s > n : Xg T ^ < r} and rj + i = inf{r > Si : X? T ^ > r}. 
Consider the killed paths := (X^^jO < t < r^ + i — Sj), i > 0, which are the excursions of 
X^ T '^ away from [r, t']. Because A^ 7- ') enters continuously in [0, r] (it has no negative jumps), 
the initial value of each of these excursions is r, except, if x < r > e o which starts from 
Because X^'^ is a killed strong Markov process, the excursions (ej)j>o form a sequence of i.i.d. 
excursions (except that eo may have a different starting point), distributed as the process X^ T > 
killed upon exiting (0, r], ending at the first excursion exiting (0, t] from the bottom. 

Then observe that A T,T = ti on [rj,Sj] and that a r,T (ij— ) = and a T,T (ti) = Sj, where 
to = and tj is the i-th hitting time of r by X( r ) defined in the theorem, so that is also the 
i-ih excursion of X^ away from r, for all i. 

As a consequence, it only remains to prove that the strong Markov process X^ T ' and the 
Levy process Y with Laplace exponent ip, both started at x < r and killed upon exiting (0, r], 
have the same law. From what precedes, we know that this law does not depend on r'. Recall 
from Theorem 13.31 that the process X^ T '^ started at x an d killed upon exiting (0, r] can be 
written 

x i t T ' ) =x-t+ Y, C(«i), o<t <t at (TiOO) , 

i:0<ti<t 

where (xj = (ui,iv(ui));i > 0) range over the leaves of C T /(T), and ti = (p{xi). By a standard 
truncation argument, we can assume that the lifespan measure A has compact support and 
choose t' sufficiently large for this support to be contained in [0,r' — r]. This guarantees that 
on [0, To A T( Tj00 )], the jumps of X^ T > are exactly the lifespans of individuals (as if no branches 
were cut down). Then the jumping rate of X( T ') from i to (x + z)dz on [0, To A T/ TjOQ \] does not 
depend on r (recall it does not depend on r' either), so we denote it by M(x,dz). Next recall 
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from the proof of the Feller property above, that the law of X^ T ) started at % and killed upon 
exiting (0, r] is the same as that of X^ T > started at x + h and killed upon exiting (h, r + h]. But 
the jumping rate at + of the former is M(x, dz), whereas that of the latter is M(x + h, dz), 
which proves that M(x, dz) does not depend on x either. As a consequence, the law of X^ T > 
killed upon exiting (0, r] is that of a Levy process (killed upon exiting (0,r]) with no negative 
jumps and finite variation, drift —1 and Levy measure, say II. It only remains to show that 
II = A. Thanks to Theorem 13.41 and Corollary 13.51 the dates at which the ancestor gives birth 
are those levels a G (0, x) such that D(a-) < D(a), where 

D(a) := inf{i > : = X ~ a G (0, X ), 

and furthermore, the lifespan of the individual born at level a is the jump size A CT of X^ T '^ 
at D(o—). Then we deduce that (A CT ;cr G (0, x)) is a Poisson point process with intensity 
measure A. But since X^ T ) is a Levy process with no negative jumps, finite variation, drift 
coefficient —1 and Levy measure II, it is known that (A a ;a G (0, x)) is a Poisson point process 
with intensity measure IT. This shows that II = A. 

(ii) Under P x (- , Ext), thanks to Theorem 14.1} we see that (i) holds for r = oo, and we 
get that under P x ( • , Ext), A is a Levy process with Laplace exponent ip, started at x, hitting 
0, and killed upon hitting 0. Now by Proposition 12.21 and Theorem 14.11 Px(^ x t) = ex P( — VX)) 
and from the Preliminaries on Levy processes, we know that this is also the probability that the 
Levy process Y with Laplace exponent A i— * A — F(X) hits 0. As a consequence, the statement 
on conditioned distributions also holds. □ 

5 New properties of splitting trees and the Crump— Mode— 
Jagers process 

In this section, we will constantly use the notation F(X) = f£°(l — e~ Xr )A(dr), and ip(X) = 
A — F(X). The scale function W is the positive function with Laplace transform (see 
Preliminaries). 

Also recall that P x denotes the law of the splitting tree with £(0) = x, whereas P x is the 
law of the spectrally positive Levy process Y with Laplace exponent ip started at x- 

5.1 Properties of splitting trees 

5.1.1 A new proof of Le Gall and Le Jan's theorem 

Theorem 5.1 Let Y be a spectrally positive Levy process with Laplace exponent ip such that 
ip(0) = and ip'(0+) > 0, started at x and killed when it hits 0. Then define 

H t = Card{0 < s < t : y s _ < inf Y r } t < T , 

s<r<t 

and L the local time process of H 

/•To 

L n := / dtl {Ht=n} n > 0. 
Jo 
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Then (L n ;n > 0) is a Jirina process with branching mechanism F starting from Lq = x> that 
is, for any integer n > 1, 

L n = S n o ••• o Sx(x), 
where the Si are i.i.d. subordinators with Laplace exponent F. 

Remark 3 In 120)] . the Levy process Y is the continuous analogue of a downwards- skip- free 
random walk that can be extracted from a Bienayme-Galton-Watson tree as follows. The jumps 
of this random walk are the offspring sizes, shifted by —1, of successive individuals of the tree 
taken in the depth-first search order. The height functional H is chosen by analogy with the 
BGW case, where it allows to recover the generation number of individuals visited in the depth- 
first search order. Then it is shown W0\ Proposition 3.2] that the time L n spent at level n indeed 
is Markovian (in n), satisfies the branching property (in its initial value), and its probability 
transitions are displayed. 

The foregoing statement is a slight refinement of Proposition 3.2 in 120], because it provides 
the law of (L n ;n > 0) as the composition of subordinators, which in passing sheds light on 
the genealogy defined in [6] by flows of subordinators. This relationship between both types of 
genealogies can be seen directly thanks to splitting trees (see proof below, and Introduction p{3|). 

Proof. Thanks to the tools set up in the previous sections, the proof of this theorem is 
straightforward. First, since ^(0) = 0, thanks to Theorem l4.1l and Proposition [221 the splitting 
tree T with lifespan measure A and law P x is subcritical (^'(0+) > 0) or critical (tp'(0+) = 0), 
and so has finite length I a.s. Now Theorem I4.3f ii) allows to state that Y is the JCCP of T, 
so in particular, To = I. 

Next, thanks to Corollary 13.51 we know that Ht is also the genealogical height of ip~ l {t) in 
T, so that 

L n = dil { |p l093 -i( t) | =n} = \({x £ T : \pi(x)\ = n}). 
Jo 

This proves that L n = J2v-\v\=nX( v ) = so the proof ends with an appeal to Theorem 14.11 
□ 

5.1.2 Exceptional points 

Recall that for any chronological tree T, 

H T = Card{x € T : P2(x) = r} < Card{x € T : P2(x) = r} = H T . 
Recall from Lemma 14.21 that for any fixed r, H r = H T < oo a.s. 
Lemma 5.2 F(Vx = (u,t) G dT, E T = oo) = 1. 

An important consequence of this lemma is that a.s. under P for all t' > r > 0, 

E T = E T = Card{t : x[ T ' ] = r}. 

In other words (H T ;0 < r < r') is a.s. equal to the occupation process of X^ T '\ Indeed, if we 
had H r < S T then there would be x £ <9T such that P2(x) = t, but in that case, S T = oo, which 
would imply H r = oo = S T . 
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Proof. Assume that there is x = (u,t) £ dT such that H T < oo. Since T has a.s. locally 
finite length, and we are interested in a local property (property of truncations), we can also 
assume, without loss of generality, that it has finite length. Recall that r = lim n uj(u\n) and 
that a(u\n) < r for all n. Since H T < oo, we must have u){u\n) < r for all but a finite num- 
ber of n's. In particular, there are infinitely many times t n := ip(u\n, uj(u\n)) increasing to 
t := tp(x) such that X tn < r = X t (and \\m. n X tn = t), where X is the JCCP of T. Setting 
U := sup{s < t : X s > X t }, we must have U < t, because otherwise, since X has no negative 
jumps, it would hit r infinitely often before time t. The bottomline is that, with positive 
probability, X has infinitely many ladder times (record times of the past supremum) on a finite 
interval. But since X is a spectrally positive Levy process with finite variation, this happens 
with zero probability [4, Chapter VII]. □ 

Actually, we have not even proved yet that dT was nonempty with positive probability under 
P. Of course, when A is finite, P(OT 7^ 0) = 0. The following statement treats the infinite case. 
Fix r > 0, and set 

r := {er € [0, r] : ~ CT = 00} and G := {s £ [0, t] : H s = 00}, 

where H s is the height of (p (s), (p^ 1 being the exploration process of C T (T), and t = A(C r (T)). 

Thanks to the paragraph preceding Theorem 13 .1\ we know that G has zero Lebesgue mea- 
sure. Thanks to Theorem 14.11 we know that T has a.s. zero Lebesgue measure. 

Theorem 5.3 If A is infinite, conditional on E T 7^ 0, T and G are a.s. everywhere dense. 

Proof. We start with G. Recall the notation X^p := X t - — Xr t _ s \, < s < t. By duality, 
and thanks to Theorem 14.3} the law of X^ on finite intervals away from and — t is that of 
the Levy process Y. Now thanks to Corollary 13.51 H s is the number of record times of the past 
supremum of X^ s \ As a consequence, to have H s = 00, it is sufficient that 

inf{r > : X^ > 0} = 0. 

But inf{r > : X^ > 0} = inf{r > s : X^ > X^}, so it suffices to show that for any finite 
interval / 

P{{s G / : inf{r > s : Y r > Y s } = 0} is dense) = 1. 

On the one hand, when A is infinite, the Hausdorff dimension of the set of times when the path 
of Y has Holder exponent h(s) < 1, is strictly positive a.s. [12]. On the other hand, it is easily 
seen that 

P(Vs G /, inf{r > s : Y r > Y s } > => h{s) = 1) = 1, 

which entails 

P(dim^({s G / : inf{r >s:Y r >Y s } = 0}) > 0) = 1, 

which in turn yields the desired result (the interval / being arbitrary). 

Now let us prove that T is everywhere dense a.s. Conditionally on H T 7^ 0, because G is 
everywhere dense and X is right-continuous, the closure of the range of G by X is [0, r]. Since 
s G G <p(s) G dT, we can conclude with Lemma 15.21 □ 
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5.1.3 The coalescent point process 

Fix t > 0. For any chronological tree T, we let (xj(r); 1 < i < H T ) denote the ranked points 
x\ <x%< ■ ■ ■ of T such that p 2 (x{) = r. In particular, the vertices pi(xi) of T are exactly the 
individuals alive at level r. 

Theorem 5.4 Conditional on H r = n > 1, let 

f p 2 (xi Ax i+ i) if ie{l,...,n-l} 
4 \ if i = n 

Then under P x (- | H T ^ 0), (aj; 1 < i < E r ) is a sequence of i.i.d. r.v. stopped at its first 0, 
whose common distribution is that of A := infYj, where Y is the Levy process with Laplace 
exponent tp started at t and killed upon exiting (0, r] . 

In particular, the duration C elapsed since coalescence between two consecutive individuals 
( if any ) has 

P(C < a) = F(A > r - a | A + 0) = \ ~_ < a < r. 

Furthermore, the coalescence level between Xj(r) and x^t) (j < k) is given by 

p 2 (xj(r) A Xfc(r)) = min{aj : j <i < k}. 
Coalescence levels can be seen on Fig. 01 
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Figure 4: Illustration of a chronological tree showing the durations ci,C2,C3 elapsed since 
coalescence for each of the three consecutive pairs (xi,X2), (x2,xs) and (x3,xa) of the H T = 4 
individuals alive at level r. 
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Remark 4 Taking A(dz) = be dz, one can recover Lemma 3 in \2J$ . Namely, since an 
elementary calculation yields W{x) = 1 + bx, 

P(Cefc) = l±^_J_ ,< T . 

For the proof of Theorem 15.41 we will need the following claim on chronological trees. 

Claim. For any chronological tree T, and any three points x, y, z E T such that x < y < z, 
then x A z ~i x A y and x A z -< y A z. In addition, 

x A z G {x A y,y A z}. (12) 

Proof. Recall from Theorem I3.4f ii) that y~<xifft<s<t. Now let x < y < z, and set 

a = <p(x), b = <p(y), c = <p(z), s = ip{x A y), and t = (p(y A z). Then since x A y ~< x,y, and 
y A z -< y, z, we get 

s < a < b < s and t < b < c < t. 

First, we show that x A z -< y. Indeed, writing r = ip(x A z), we have r < a < c < r, but since 
a < b < c, we get f < b < r, which is exactly x A z -< y. Then notice that since x A z ~< y and 
x A z -< x, z, we always have x A z ~< y A z and x A z ~< x A y. This is the first assertion of the 
Claim. 

Next, we show that ([I2|) <^> (Tl3]) . where 

x Ay < z or y A z ^ x. (13) 

Assume (fT3l) . If x A y -< z, then since x Ay ^ x, we get xAj/ ^ iAz; on the other hand, thanks 
to the first assertion of the Claim, x A z -< x Ay, so that x A z = x Ay. The same reasoning 
shows that if y A z -< x, then x A z = y A z. Thus we have shown (fT3|) => (JT2J) , and the converse 
implication is straightforward. 

Now observe that (fT3l) is equivalent to 

s<c<s or t < a < t. 

To prove (|12p . we assume that the last assertion does not hold, and show that this will bring 
a contradiction. Indeed, if the last display does not hold, then c > s and a > t, so that 

s<a<i<b<s<c<t. 

In particular, we can extract from the last display the following two double inequalities s <t < s 
as well as t < s < t, which reads also 

ip-\t) <^-\s) <^- l (t). 

Writing ip~ 1 (s) = (u,o~) and 99" 1 (t) = (v,t), and recalling from Theorem I3.4f i) that ip~ 1 (s) = 
(u,Ld(u)) and y>~ l {t) = (v,uj(v)), we see that (u,a) is in the segment [(v,t), (v,u(v))], from 
which we conclude that u = v. Therefore, s = t, which contradicts the fact that s < a < t. □ 
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Proof of Theorem 15.41 . Let us start with the last assertion of the theorem. A straight- 
forward consequence of the Claim is that for any x < y < z, we have pzix A z) = mhx(p2{x A 
y)>P2(y A z)). Indeed, from x A z -< x A y and x A z -< y A z, we get Az) < min(p 2 (^ A 
y)iP2(y A z)), and the converse inequality stems from (j 12 [> . The last assertion of the theorem 
then follows from a recursive application of this property. 

Next, recall from Theorem 14.31 the following recursive definition: to = and = inf{t > 
ti : G {0, t}}. Then abusing notation (confounding those for T and C T (T)), we have 

ti = <p(xi), and by Theorem I3,4f iii) . we get that for any 1 < i < E T , <Zj is given by 

at = inf X t (r) . 
ti<t<ti_)_i 

But thanks to Theorem l4.3l under P x (- | S r 7^ 0), the killed paths ej := (A^ 1 ^, < i < ti+%— ti), 
i > 1, form a sequence of i.i.d. excursions, distributed as the Levy process Y started at r, 
killed upon exiting (0, r], ending at the first excursion hitting before (r, +00). It is then 
straightforward that (dj; 1 < i < E T ) is a sequence of i.i.d. r.v., stopped at its first 0, and 
distributed as A, where, thanks to the Preliminaries on Levy processes, 

F(A < a) = P T (Y exits (a,r] at the bottom) = W(0)/W(t - a) = 1/W(r - o") cr G [0,r]. 

The fact that W(0) = 1 can be deduced from a Tauberian theorem, checking that 
converges to 1 as A — > 00. In particular, the duration C = r — A elapsed since coalescence 
between two consecutive individuals has 

rfC <a) nA>r-a) _ l- 1/W(a) 

n °- j " P(^0) "1-1/W(r)' 

which completes the proof. □ 



5.2 Ages and residual lifetimes 

We hope that the last subsection has convinced the reader that multiple uses of Theorem 14.31 
can be made by applying standard results on Levy processes to one's favourite question of 
population biology. Hereafter, we give a simple example, in the form of a statement on ages 
and residual lifetimes of living individuals. The proof of this statement is left to the reader, 
since it relies on the same analysis as done in that of Theorem 15.41 

Fix r > 0. As in the previous subsection, we let (xj(r);l < i < 3 T ) denote the ranked 
points x\ < x% < ■ ■ ■ of T such that P2(xi) = r. Set also Uj := pi(xi) the z-th individual of T 
who is alive at level r. Then for 1 < i < H T , 

Ai := t — a(ui) and Ri := u>(ui) — r 

define respectively the age and residual lifetime of Uj at 'time' r. 

Proposition 5.5 Under P x , conditional on H T = n > 2, the individuals {ui\2 < i < n) alive 
at time r except u\, have i.i.d. ages and residual lifetimes, whose common distribution is 
independent of x an d equal to that of a pair (A T ,R T ) such that 

F(A T e dx, R T G dy) = Pq{-Y T{0 oc] _ G dx, Y T{0 oo] G dy \ T (0;Oo] < T_ r ) < x < r, y > 0, 

where Y is the spectrally positive Levy process with Laplace exponent i/j. 
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5.3 Properties of the Crump— Mode— J agers process 
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Figure 5: Graph of the Laplace exponent ip of the Levy process Y (in the finite case). On this 
example is shown a supercritical exponent, with positive largest root n, such that ip(0) = 
(i.e. A({+oo}) = 0). The root rj is the Malthusian parameter of the CMJ process. In the finite 
case, note that because tp is convex, and tp(\) > A — b, one has r] < b. 

Proposition 5.6 The probability of extinction is P x (Ext) = e _,?x . In addition, the one- 
dimensional marginal of S is given by 

¥ x (E t = 0) = W(t- X )/W(t), 

and conditional on being TiOTizerOj l^j T has a geometric distribution with success probability 
1/W(r). In particular, E x (~ r | E T ^ 0) = W(t). 

Proof. The value of the extinction probability stems straighforwardly from Proposition 12.21 
and Theorem 14.11 

Since H r = iff the first excursion of X^ T ' hits before (r, oo), Theorem I4.3( i) entails 

F X (E T = 0) = P X (T < r (T>+oo) ) = W(t - x)/W(r). 

Now thanks to Theorem I4.3l fi) again, conditional on being nonzero, H T is the length of a 
sequence of i.i.d. excursions of Y starting from r and killed upon exiting (0, r], stopped at the 
first one that exits it from the bottom. As a consequence, the conditional distribution of H r is 
geometric with success probability 

P t (To<T (Ti+oo) ) = 1/W(t). 
The computation of E x (H r | E T ^ 0) is then straightforward. □ 

Proposition 5.7 Set P* 1 := P(- | Ext). Underf \ the (supercritical) splitting tree has the same 
law as the subcritical splitting tree with lifespan measure e~ r,r A(dr) . In particular, in the finite 
case, its birth rate equals b — n. 
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Proof. Thanks to Theorem I4.3f ii). we know that under F\ the JCCP X is a spectrally 
positive Levy process with Laplace exponent ip conditioned on hitting 0. Now it is known [1] 
that this conditioned Levy process is a spectrally positive Levy process with Laplace exponent 
ip\ where tfr{\) = + tj). It is straightforward to obtain the following equality 

ifl(\) = \- e'^il - e~ Xr )A(dr) A > 0, 
Jo 

which proves that under P^, X is the JCCP of a splitting tree with lifespan measure e~ r,r A(dr). 
It is indeed a subcritical splitting tree, since ^(0) = ip(>q) = and ^'(0+) = if/(rj) > 0. 
Another elementary calculation shows that in the finite case, the birth rate $ = J °° e - ^ A(dr) 
of this splitting tree equals b — r\. □ 

Proposition 5.8 The asymptotic behaviour of 3 is as follows. 

(i) (Yaglom's distribution) In the subcritical case, 

lim P(~ r = n | H T ^ 0) = m n_1 (l - m) n > 1. 

T—>00 

(ii) In the critical case, provided that J°° r 2 A(dr) < oo, 

lim P(H r /r > x \ E T / 0) = exp(-/'(0+) x/2) x > 0. 

r^oo 

fmj In i/ie supercritical case, conditional on Ext c , 

lim e^ T E T = f , 

where £ is an exponential variable with parameter ip'{rj). 

Remark 5 Whereas the convergence stated in (Hi) certainly holds solely in distribution when 
A is infinite (because of Theorem \5.3\) . there is actually pathwise convergence in the finite case, 
provided that J°° r log(r)A(dr) < oo [22]. Thanks to this result, rj is seen to be the so-called 
Malthusian parameter. 



Proof. We know from Proposition 15.61 that conditional on being nonzero, H r is geometric 
with success probability 1/W(r). 

In the subcritical case (ip(0) = and ip'(0 + ) = 1 — m > 0), ip(X)/A converges to 1 — m as 
A vanishes, so a Tauberian theorem entails that W(t) converges to 1/(1 — m) asr-> oo. This 
proves (i). 

In the critical case (ip(0) = and tp'(0 + ) = 0), ^(A)/A 2 converges to t/>"(0 + )/2 as A vanishes, 
so a Tauberian theorem entails that W(t)/t converges to 2/ip"(0 + ) as r — > oo. This proves 
(ii)- 

As for (iii), set W^x) := W(x)e r)X . Then the Laplace transform of is 

/•OO 1 

/ W\x)e- Xx dx = — A > 0, 

Jo ^(A + r?) 

which is equivalent to 1/^'(t])\ as A vanishes. As previously, we deduce that W\t) converges 
to l/il)'(rj) as r — > oo, which reads 

lim W{r)e-^ = 
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The convergence in distribution then follows from the asymptotic equality between P(H r e VT € 
• | E T ^ 0) and P(H T e" , ' T € • | Ext c ). □ 

Actually, we can go a little further than the previous proposition, by distinguishing between 
points whose descendance is either finite or infinite. Recall that finiteness here refers to the 
first projection of the tree (discrete part), and set 

3^° := Card{x G T : pz(x) = r,8(x) is infinite}, 

3f := Card{x E T : P2(x) = t, 6{x) is finite}. 

In particular, + E£ = S r . 

Proposition 5.9 Set p := ^'{rf) < 1. In the supercritical case, conditional on Ext c , 

lim e-v (S~s£) = (pt,(l - p)H), 

r^oo 

where £ is an exponential variable with parameter p. In particular, e~ ril "E'^ converges in dis- 
tribution to an exponential variable with parameter 1. 

Remark 6 In a work in preparation \1T$ , we show that actually, (H^°;r > 0) is Markovian, 
and, more precisely, that it is a Yule process with birth rate r\. In particular, the convergence 
of e~ VT to an exponential variable with parameter 1 is known to be pathwise. 

Proof. We use the notation of Proposition 15.51 Each individual U{ alive at level r has infinite 
descendance iff the splitting tree starting from m x (r, u(ui)) is infinite, which, conditionally 
on the residual lifetime Ri := u)[ui) — r, occurs with probability 1 — e~ vRi . Now thanks 
to Proposition 15.51 conditionally on H T = n > 2, the residual lifetimes of the individuals 
{ui\2 < i < n) are i.i.d., all distributed as R T . As a consequence, by the branching property, 
the conditional distribution of (2j?°,Sr), modulo the first individual u\ (which is negligible as 
the population size goes to infinity), is that of (B,n — 1 — B), where B is a binomial variable 
with parameters p T and n — 1, and p T := E(l — e" 71 ^). Conditional on non-extinction, the 
only change is that B is conditioned to be greater than 1. Now since the population size goes 
to infinity as r — > oo, this last conditioning vanishes, but above all, (H^° , s£)/H T converges 
in probability to (poo>l — Poo), with p^ := lim r p r . Therefore, it only remains to show that 
Poo = P, since we know from Proposition 15.81 that e~ r,T 'B T converges in law to £. Now recall 
that p T = E(l — e~ vRT ), and observe that R T converges in distribution, as r — > oo, to the law 
of Yt (0 under Pq. Summing up, we have to prove that 

p = E (l-exp(-r)Y T(0oo] )y 

As a first step, we apply the compensation formula to the Poisson point process (A s ; s > 0) of 
jumps of Y. Denoting by St the past supremum of Y at time i, we get, for any A > 0, 

E (exp ( - AY T(0)Oo] )) = E 1 {s t .=o} exp ( - A(Y t _ + A<))l { y t _ +At>0} 

t:At>0 

roo poo 

= E dt A(dr) l {St =o} exp ( - X(Y t + r))l {Yt+r>0} 
Jo Jo 

roo rr 

= / A{dr)e~ Xr / fi Q (dx)e Xx , 
Jo Jo 
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T(0,oo] 



where we have set 

Ho(dx) := E I dtl { _ Yt edx} x > 0. 
Jo 

In order to compute an expression for [j,q, we also define 

H(dx) := E / dt l{_ Yt adx} x > 0. 
Jo 

Because Y has derivative a.s. equal to —1 when it hits — x, we deduce that fi(dx) = Eq(N x ) dx, 
where N x is the number of passage times at —x. As a consequence, applying the strong Markov 
property at T_ x , we get (see Preliminaries) 

fi(dx) = P (T_ X < oo)E_ x (N x ) dx = E {No)e- vx dx x > 0. 

On the other hand, applying the strong Markov property at the successive hitting times of 
under Pq, say = Si < S2 ■ ■ ■ < Sn < Sn +i = 00, we get 

N o rS n+1 

H(dx) =E ^2 dt 1 {-Y t edx} = E (N )fi (dx), 

n=l J Sn 

so that, in conclusion, 

Ho(dx) = e _7?a: dx x > 0. 

Therefore, for any A > 0, 

E (exp ( - XY T(0 J) = ^ A(dr)e- Xr jf e^>, 

so in particular, 

E (exp ( - TjV Tl0iOo] )) = A(dr) r e^ r = 1 - ^(tj), 
which is the desired result. □ 

The Markovian (finite) case. There are two possibilities for (H T ;r > 0) to be Markovian. 

First, when A is a Dirac mass at {00} (wih mass 6), (S r ; r > 0) is a pure-birth process (with 
rate b), and since tp(X) = X — b, we get W(x) = exp(6x), so that the size of the population at 
'time' t is (shifted) geometric with success probability exp(— br). 

Second, when A is exponential, with parameter, say d, then (H r ;r > 0) is a birth-death 
process with birth rate b and death rate d (supercritical iff b > d), so that A(dr) = bde~ dr dr. 
Then we have 

In particular, rj = b — d and ip'(r]) = 1 — (d/b). The scale function (see also Remark |4]) is 
W[x) = 1 + bx in the critical case (d = b), and in all other cases, 

TTr/ d-be^-V* 

W(x) = — x > 0. 

One recovers that conditional on being nonzero, the size of the population at 'time' r is geo- 
metric with success probability (d — b)/{d — be^ b ~ d ^ T ). In the supercritical case, conditional on 
non-extinction, the fraction of individuals with finite descendance converges in probability to 
d/b, whereas the number of individuals with infinite descendance is a Yule process with birth 
rate b — d. 
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